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Abstract 

By generalizing the Green's function approach proposed by Beliaev Olli we investigate the 
effect of quantum depletion on the energy spectra of elementary excitations in an F = 1 spinor 
Bose-Einstein condensate, in particular, of 87 Rb atoms in an external magnetic field. We find 
that quantum depletion increases the effective mass of magnons in the spin-wave excitations 
with quadratic dispersion relations. The enhancement factor turns out to be the same for both 
ferromagnetic and polar phases, and also independent of the magnitude of the external magnetic 
field. The lifetime of these magnons in a 87 Rb spinor BEC is shown to be much longer than that 
of phonons. We propose an experimental setup to measure the effective mass of these magnons 
in a spinor Bose gas by exploiting the effect of a nonlinear dispersion relation on the spatial 
expansion of a wave packet of transverse magnetization. This type of measurement has practical 
applications, for example, in precision magnetometry. 

Keywords: Spinor Bose-Einstein condensates (BECs), Beliaev theory, Energy spectrum, Spin 
wave, Beliaev damping 



1. Introduction 

Since the experimental realization of Bose-Einstein condensates (BECs) iS 3. IH, the Bo- 
goliubov theory of weakly interacting dilute Bose gases has been successfully applied to de- 
scribe a variety of phenomena in these systems flfjl LZl |8|. The Bogoliubov theory was originally 
invented to describe bosonic systems at absolute zero H, and then extended to finite tempera- 
ture It gives the leading-order values of physical observables of a system in 
thermodynamic equilibrium. The second-order correction to the Bogoliubov result is usually 
relatively small for weakly interacting dilute Bose gases. At absolute zero, this correction is a 



consequence of a small fraction of quantum depleted noncondensed atoms 114111511 . The second- 
order correction to the Bogoliubov energy spectrum was given by Beliaev who developed a 
diagrammatic Green's function approach to describe the energy spectrum of elementary excita- 
tions at absolute zero Afterwards, finite-temperature theories based on the Beliaev technique 
were developed for weakly interacting Bose gases lflol[l6l 17, 18]. With the rapid development 



of techniques for precise measurements of physical observables, the small effect of quantum 
depletion is no longer beyond the scope of experimenters. Furthermore, by using a Feshbach 
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resonance or optical lattices, the effective interatomic interaction can be manipulated to cover 
both weakly and strongly interacting systems £01 III HQ. In particular, it has been shown 
that up to a moderate strength of interaction, by taking the second-order correction to the mean- 
field (Bogoliubov) calculation, the obtained results for spinless condensates agree excellently 
with both the results of experiment and those of quantum Monte-Carlo simulation 112211 . There- 
fore, the second-order correction to the Bogoliubov result, which can be obtained in an analytic 
form, can be used as an important check for any calculation or measurement of a strongly corre- 
lated system. Furthermore, the Beliaev theory also predicts the so-called Beliaev damping which 
quantitatively shows a finite lifetime of Bogoliubov quasiparticles (phonons) due to their colli- 
sions with condensed particles. The Beliaev d amping of quasi-particles under various conditions 
has been a subject of active study 11231 l24i l25l \26l 12711 . 

Recently, Bose-Einstein condensates with spin degrees of freedom (spinor BECs) have been 
extensively studied (see, for example, ll28lp . These atomic systems simultaneously exhibit su- 
perfluidity and magnetism, and the combination of atoms' motional and spin degrees of freedom 
gives rise to various interesting phenomena in the study of thermodynamic properties and quan- 
tum dynamics. Due to the competition between spin-dependent interatomic interactions and the 
coupling of atoms to an external magnetic field, the system can exist in various quantum phases 
with different spinor order parameters [29 30j 31 1. In contrast to spinless BECs, there exist 
spin-wave excitations in spinor BECs in addition to the conventional density-wave excitations. 
These are excitations of atoms from the condensate to the other magnetic sublevels, and the cor- 
responding magnons have quadratic dispersion relations at low momenta as opposed to the linear 
dispersion relations of phonons. Furthermore, in spinor Bose gases the collisions of atoms in dif- 
ferent spin channels give rise to spin-conserving and spin-exchange interactions. Particularly, 
in some atomic species such as 87 Rb, the ratio of the spin-conserving to spin-exchange interac- 
tions is so large that it can compensate for the small noncondensate fraction. That is, the mean 
field caused by noncondensed atoms with spin-conserving interaction can have the same order of 
magnitude as that caused by condensed atoms with spin-exchange interaction. Consequently, a 
small number of noncondensed atoms can, in principle, give an appreciable effect on the physical 
properties of the system, for example, by shifting the phase boundary between different quantum 
phases 113211 . 

In this study, we apply the Beliaev theory to spin-1 Bose gases to investigate the effect of 
quantum depletion at absolute zero on the energy spectra of elementary excitations. In the pres- 
ence of an external magnetic field, the ground state can be in several quantum phases, depending 
on the strength of the quadratic Zeeman energy relative to the spin-exchange interatomic inter- 
action. In contrast to the work in |32|], we do not consider phase transitions between different 
quantum phases. Instead, we assume that the magnitude of the external magnetic field is chosen 
so that the system is stable in a certain quantum phase. Here, we consider two characteristic 
phases of F = 1 spinor Bose gases: the fully spin-polarized ferromagnetic phase and the unmag- 
netized polar phase. In the calculation of second-order corrections for ultracold atomic systems 
like 87 Rb, the spin-conserving interaction must be taken into account while the spin-exchange 
interaction is neglected because of its much smaller value. (The spin-exchange interaction is, of 
course, taken into account in the calculation of first-order values.) We find that for both the fer- 
romagnetic and polar phases, the quantum depletion leads to an increase in the effective mass of 
magnons, while it does not alter the energy gap to the leading order. Although the effective mass 
is different between the ferromagnetic and polar phases, it is enhanced by the same factor for 
these quantum phases. This factor is also independent of the magnitude of the external magnetic 
field. This implies a physical mechanism whereby the quantum depletion affects the motion of 
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quasiparticles in spinor Bose gases in a universal manner under some certain conditions. In the 
case of 87 Rb, where the spin-conserving interaction is much larger than the spin-exchange one, 
the lifetime of magnons becomes much longer than that of phonons. We show that this agrees 
with the mechanism of Beliaev damping which is caused by collisions between quasiparticles 
and the condensate. To measure the effective mass of magnons in spinor Bose gases, we propose 
an experimental scheme which exploits the effect of a nonlinear dispersion relation on the spatial 
expansion of a spinor wave packet during its time evolution. This type of measurement can be 
used for several applications: to probe the effect of quantum depletion, to identify spinor quan- 
tum phases, or to be used for precision magnetometry in a way different from the method given 



This paper is organized as follows: Section [2]formulates the diagrammatic Green's function 
approach for spin-1 spinor BECs, which is the generalization of Beliaev theory to systems with 
spin degrees of freedom. The explicit forms of the matrices of self-energies for both the ferro- 
magnetic and polar phases are given. The T-matrix which plays the role of an effective interaction 
potential in dilute Bose gases is also introduced in this section. Section[3]summarizes the results 
of energy spectra of elementary excitations at the first order in the interaction. It is the rederiva- 
tion of the Bogoliubov energy spectra by using the Green's function approach ll34ll . Section [4] 
deals with the self-energies to the second order in the interaction and gives the leading-order 
corrections to the Bogoliubov energy spectra due to the effect of quantum depletion. Section [5] 
shows that the elementary excitations with quadratic dispersion relations are spin waves. An 
experimental scheme using spinor wave packets is proposed to measure the effective mass of 
magnons. An order-of-magnitude estimation of the time evolution of these wave packets is also 
given in this section. Section[6]concludes the paper by discussing the application of the measure- 
ment to some practical purposes. The detailed calculations are given in the Appendices to avoid 
digressing from the main subject. 

2. Green's function formalism for a spinor Bose-Einstein condensate 

2.1. Hamiltonian 

We consider a homogeneous system of identical bosons with mass M in the F = 1 hyperfine 
spin manifold that is subject to a magnetic field in the z-direction. The single-particle part of the 
Hamiltonian is given in the form of a matrix by 



where the subscripts j,f — 0, +1 refer to the magnetic sublevels, and cjb is the coefficient of 
the quadratic Zeeman energy. Because of the conservation of the system's total longitudinal 
magnetization, the linear Zeeman term vanishes. The total Hamiltonian of the F = 1 spinor Bose 
gas is then given in the second-quantized form by 



in [|33ll. 



h 2 v 2 



+ lBj 2 8jf, 



(ho)j f = - 



2M 



(1) 




if 



(2) 



where 0j(r) is the field operator that annihilates an atom in magnetic sub-level j at position r, 
and the interaction energy "V is given by 




j,j',m,m' 



(3) 
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Here, the matrix element V} m j/ m <(r - r') can be written as a sum of interactions in two spin 
channels T — and 2 (T denotes the total spin of two colliding atoms) as follows: 

V jm , r „Ar - r') = (j, m\T = 0)(T = 0\f, m')V (r - r') 

+ (j, m\r = 2){T = 2\f, m')V 2 (r - r'), (4) 

where quantum statistics prohibits bosons from interacting via the spin channel T = 1. 

In the presence of a condensate, the field operator 0j(r) is decomposed into the condensate 
part, which can be replaced by a classical field V«o^> and the noncondensate part 6j(r): 

f,(r) = + Sj(r). (5) 

For a homogeneous system, the condensate is characterized by the condensate number density 
no and the spinor order parameter gj(J = 0, + 1 ), which is normalized to unity: 

!>/ = !• (6) 
Substituting Eq. (01 into Eq. (0, we can decompose the interaction energy as 

7 

<V = E + V„, (7) 



where 



(8a) 
(8b) 



Eo=\n 2 J dr J dr't*CVj,n,r,Ar-r% n ,{ r , 
t>i =\n j dr J drTjCVjmj> m '(r-r')6 m ,(r')6 f (r), 
%=\n G J dr J dr'5]:(r)^(r')V; mj w(r-r')^^, (8c) 
% =2(i« ) J dr J" dr'^r')V M ,w(r - r')^^(r), (8d) 
t>4=2(i« ) J dr J dr'«j(r)^y >li/j/ (r-r'fc* / (r), (8e) 
V 5 =2(^n l l 2 ) jdrj dr'^W^CrOV^jwCr-r'^^Cr), 
V 6 =2(|«y 2 ) J"dr J dr'^(r)^y M/m ,(r-r')^(r')^(r), 
t/ 7 =i J dr J" dr'^(r)^(r')F M/m '(r-r')^(r')§/(r). 



(8f) 
(8g) 
(8h) 



These interactions are illustrated by the Feynman diagrams in Fig. Q] 

We consider a grand canonical ensemble of the above atomic system, and introduce the op- 
erator 

K^ii-liN, (9) 
4 



N J 

V 



m / 



NAAMA/V 



' J" 




m 



m 




f 



m 





A 

v 7 



Figure 1: Two-particle interactions involving different numbers of condensed and noncondensed atoms. The dashed, 
solid, and wavy lines represent a condensed atom, a noncondensed atom, and the interaction, respectively. 



where ju denotes the chemical potential and TV is the total number operator: 

/V= f dr^J(r)^(r). 

J j 

Using Eqs. (|2]i,(|5]),(|7]), and (0, we have 

r 



(10) 



7C = £o 



No+K', 



(11) 



where Eq, given in Eq. ( I8al ). is the interaction energy between condensed atoms, A^o = Vno is the 
total number of condensed atoms with V being the volume of the system, and 



is the corresponding operator for the noncondensate part with 

%) = (4 -V + qBj 2 )a) M ajx, 
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(12) 

(13) 
(14) 



Here, e£ = /i 2 k 2 /(2M) is the kinetic energy of a particle with momentum fik, and a v j< is related 
to the noncondensate field operator o/(r) via a Fourier transform: 



(15) 



In the following sections, "TCo and < K\ are referred to as the noninteracting and interacting parts 
of operator TC in Eq. JT21 . respectively. For a weakly interacting system, "7Ci can be treated as a 
perturbation to < K§. 



2.2. Green 's functions 

In the presence of the condensate, the Green's function is given by |2, 351 



iG t ?p x (x,y) = no^e f + iG jf (x,y), 



(16) 



where j, f — 0, +1 indicate the spin components, and x = (r, t),y = (r', f') are four-vectors in the 
time-coordinate space. The noncondensate part of the Green's function is defined as 



iGjy{x,y) = 



(O|'T^hW^,h0 ; )I O > 



(17) 



Here, |0) is the ground state of the interacting system, and T and H denote the time ordering 
operator and the Heisenberg representation, respectively. 

In the presence of the condensate, we must take into account the collision processes in which 
two noncondensed atoms get into or out of the condensate. For this purpose, in addition to the 
normal Green's functions Gjy{x,y) defined in Eq. ( TTTT i. it is necessary to introduce the so-called 
anomalous Green's functions which are defined as 



(Q\T6] H (x)6], H (y)\0) 

= jo\o } . 

. (OIT-^hM^hOOIO) 
= <5jO> ' 



(18) 
(19) 



In energy-momentum space, the Dyson's equations for the noncondensate Green's functions 
are given by 



GpP) = (G%(p) + (GX^AP)G% r (p), 



(20) 



where hp = h(po, p) is the four-momentum, and a,/3, y, 6 — 1,2 are used to label the normal and 
anomalous Green's functions as matrix elements of a 6 x 6 matrix: 



G\\(p) G\\{p) G\\ X { P ) G|» G\ 2 Q (p) 
G l Q \( P ) Gj» G'^O?) G'» G l Q 2 Q ( P ) 



G^ip) 



G l _\ x (p) G u l0 (p) G'X^p) G n hl { P ) G%(p) G'X^ip) 



G\\( P ) Gf () ( P ) Gf_,(p) Gf» Gf (p) 
Gl\(p) G 2 l (p) G^ip) Gl\(p) G 22 (p) 
G 2 \,{p) G 2 \ (p) G 2 \ _j(p) G 22 x l (p) G 22 ( P ) G 22 ^) 



G 2 1,(P) 



(21) 



where 



G)),(p) = G jr (p), G 22 ( P ) = G, f {-p). 



(22) 



Equation (1201 . which is illustrated in Fig. [2] can be written in terms of 6 x 6 matrices as a 
matrix equation: 



G(p) = G\p) + G°(p)±(p)G(p), 



(23) 




Figure 2: Dyson's equations for the normal and anomalous Green's functions. The thick line, thin line, and oval represent 
the interacting, non-interacting Green's functions, and the proper self-energies, respectively. 
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where G, G°, and £ denote the 6x6 matrices of Green's functions, non-interacting Green's 
functions, and proper self-energies, respectively. The normal and anomalous self-energies are 
labeled in the same way as the Green's functions. The solution to Eq. ( 1231 ) can be written formally 

as 



G(p) = [l-G (p)±(p)]~~ 1 G°(p). 
The non-interacting Green's function is defined as 

<o|rS;,H (*)SJ, iHo Cy)|o> 



<0|0> 



(24) 



(25) 



where |0) is the non-interacting ground state, and Ho indicates the free time evolution in the 
Heisenberg representation under the non-interacting Hamiltonian TCq given by Eq. dI5V Here, 
|0) is the vacuum state with respect to noncondensate operators; that is, flk,;|0) = for all k + 
and j = +1,0. Substituting Eq. (fT3l l into Eq. d25l l. we obtain the Fourier transform of Gj.,(x - y) 
as 



G>) = JV*<T^G° /( *) 



■■Si, 



p {) - e°/h + fi/H - q B j 2 /h + ir\ 



= S Jf Cfj(p), 



(26) 



where rj is an infinitesimal positive number. Note that the anomalous Green's functions in a non- 
interacting system are always zero, and thus, the matrix G a (p) is diagonal with matrix elements 
given by Eq. ( 1261 ). 

Now we consider two cases in which the mean-field ground state is in the ferromagnetic 
phase and in the polar phase. 

2.2.1. Ferromagnetic phase 

If the system's ground state is in the ferromagnetic phase, the condensate's spinor is given by 



(ft,6,f-i) = (1,0,0); 



(27) 



i.e., all condensed atoms reside in the 7 = 1 magnetic sublevel. Then, the only nonzero matrix 
elements of t,(p) are 



















2J!o(P) 


















































^o,o( _ 


























(28) 



This can be understood by considering the spin conservation in normal and anomalous self- 
energies, which are illustrated in Fig. [3] For normal self-energies Z'.J,(p)> the conservation of 
the total projected spin allows only j = /, i.e., diagonal elements. In contrast, for anomalous 




Figure 3: Normal and anomalous proper self-energies of noncondensed particles, hp is the four-momentum, while j, f 
label the spin components. The dashed lines represent condensed particles. 



self-energies 2,\y(p), only the j = f = 1 element is nonvanishing because the condensed atoms 
are all in the mp = 1 magnetic sublevel. 

By substituting Eq. d28l i into Eq. (f24t and using the fact that G°(p) is a diagonal matrix (see 
Eq. d26b~). we find that the matrix G(p) of interacting Green's functions has the same form as 
t(p): 

\G u (p) G^ip) 

Go.oO) 

G_!_i(p) 

G]\(p) G U H>) 

Go fi (-p) 

G_i,_,(~p)J 



(29) 



Both G(p) and £(/;>) are block-diagonal matrices composed of one 2x2 and four lxl sub- 
matrices. 

The normal and anomalous Green's functions given by Eq. (124} then can be expressed in 
terms of the self-energies as 



G U (P) = 
Go,o(P) = 



-[GO(-p)]- 1 + l\\(-p) po + e a p /H + q B /h + X\\(-p)-fx/h 



1 



[G»]- I -S'» + /^ 



, G-i-i(p) = 



[G° 1 (p)]- 1 -Ei 1 1 _ 1 (p) + i» 7 ' 



D, 



Gfi(p) 



(30a) 
(30b) 

(30c) 
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where 



Di = - [Glip^lG^i-pT 1 +X\\(p)[G° 1 (-p)r 1 +X\\(-p)[G° l (p)T 1 
- Z\\(p)Z\\(-p) + X^ip^ip) + ir, 

= pl- - 1\\{- P )\ po + *\\{p)^{p) 

-[e p °/ri-W^ + W^+^ ^ ] +(— 2^ ) (3D 

From Eqs. ( f30b and (l3TT l, we obtain the modified version of the Hugenholtz-Pines condi- 
tion J36ll for an F — 1 spinor BEC in the ferromagnetic phase, that is, for the three elementary 
excitations to be gapless, the following condition must be met: 



E)> = 0, p = 0) - E]> = 0,p = 0) = (41 - q B f)/H. 



(32) 



Here, the excitation modes with spin j = 0, - 1 are single-particle like, and thus, the correspond- 
ing anomalous self-energies and Green's functions vanish. The energy shift of -qg from the 
chemical potential on the right-hand side of Eq. d32l results from the difference in quadratic Zee- 
man energy between magnetic sublevels j = ±1 and j = 13711 . For the ferromagnetic phase, the 
Hugenholtz-Pines condition ( 1321 holds only for j = 1 in the presence of the quadratic Zeeman 
effect; therefore, only the corresponding phonon mode (j = 1) is gapless. When q^ = 0, the 
spin-wave mode ( j = 0) also becomes gapless with a quadratic dispersion relation. 



2.2.2. Polar phase 

If the system's ground state is in the polar phase, the condensate's spinor is given by 

£-1) = (0,1,0); 



(33) 



that is, all condensed atoms occupy the j = magnetic sublevel. With an argument similar to the 
ferromagnetic phase , the only nonzero matrix elements of £(p) and G(p) are the following: 








G u {p) 







-1,1 



ip) 











Go,o(p) 



Gf fi (p) 







£ n i,-i0>) 







G-iMP) 
Gf_ x {p) 










sftc-lo 







Gnl-P) 







^ao(p) 







G^ip) 



Go,o(~p) 




-10>) 











G-u-d-p) 



(34) 



(35) 



Both of these matrices are block-diagonal matrices composed of three 2x2 sub-matrices. Here, 
Sj 2 j(p) andGj 2 j(p) are nonzero due to the projected-spin-conserved scattering process in which 
two condensed atoms both in the spin state 7 = collide with each other to produce two noncon- 
densed atoms with spin components j = ±1 (see Fig. [3}. 
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The normal and anomalous Green's functions given by Eq. (l24t can then be expressed in 
terms of the self-energies as 



G u (p) = 




/ \ 0/7- / 7- X"" 1 1 1 

-i(-.P) Po + £p/n. + q B /n + l.i\ _, 


(-p)-n/h 


(36a) 








Go.oO) = 




-p) p + 4/H + ^ (-p)-fi/h 




(36b) 


D 


Do 




G-i,-i(p) = 


-[Cfi 1 (.-P)r 1 +T,\\(- 


-P) po + e^/h + q B /fi + I.\ l l (- P ) 


-p/h 


(36c) 


D-i 


D-i 




g\Up) = 




-^?.-iCp) 




(36d) 


g\Up) = 


Do ' G o> )= ' 






(36e) 


G-Up) = 


- Z -l 2, 

D i > G -> ): 






(36f) 



where 

Oi =-[G?(p)]- 1 [G° 1 (-p)]" 1 +S} j 1 1 (p)[G° 1 (-p)]- 1 +Si 1 1 _ 1 (-p)[G?(p)]- 1 

- z^p^^-p) + E^oozl^oo + 6? 

= p2 _ ^n (p) _ ^(-p)] jx, + 22i,i(P)Slf_i(p) 

" [ S + ^ 1 + ( 2 ) + *> < 37a > 

D Q =- [Gg(p)]- 1 [Gg(-p)]- 1 +Z» 1 (p)[G°(-p)]- 1 +4 1 (-p)[G°(p)]- 1 

- SjJ,(p)ZJi(-p) + I» (p)Sg,(p) + »7 

= p 2 - [^(p) - 4!o(-P)] Po + ^ (pKo(p) 

" [— + 2 1 + ( 2 ) + (3?b) 

D-i =- [G° 1 (p)]- 1 [G?(-p)]" 1 +2[j'- 1 (p)[G?(-p)]- 1 +2| 1 1 (-p)[G° 1 (p)]- 1 

- £ii,_ 1 (p)Zl > 1 1 (-p) + Z? 1 _ 1 (p)Si 2 11 (p) + //? 

_ _ Sj\(-p)]p + ^(p^iip) 

t 4-f + <lB gU(p)+^(-p) 2 / Si 1 1 ,_ 1 (p)-4 1 1 (-p) 2 . 

- [ s + 2 ] + ( ) + m. (37c) 

From Eqs. (l36l l and (l37l i. we obtain the modified version of the Hugenholtz-Pines condition 
for an F — 1 spinor BEC in the polar phase, that is, for the three elementary excitations to be 
gapless, the following condition must be met: 

Ej,' (po = 0, p = 0) - Zf_/p = 0, p = 0) =(p - qshlh. (38) 

For the polar phase, the Hugenholtz-Pines condition (l38l holds only for j — in the presence of 
the quadratic Zeeman effect (qg + 0); therefore, only the corresponding phonon mode ( j — 0) is 
gapless. 
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2.3. T-matrix 

For a weakly interacting dilute Bose gas, the contributions from all ladder-type diagrams to 
the self-energies are shown to be of the same order of magnitude [[l],[35l]> ar >d, therefore, all of 
these contributions must be taken into account. The T-matrix is defined as the sum of an infinite 
number of ladder- type diagrams as illustrated in Fig. |4] It is written as 

^jmJ'm'iPUPi; Pi, PA) - Vjmj'm'ipi - P3) 

j",m" " 

x V;,„,/<m"(q)V rm », /m <(pi - q - P3) 
+ . . . 

= V w (Pi-P3)+ 2, 7^)3 

j".m" J 

1 



Hpi)o + HP2)o - 4,- q - fp 2+q + 2/i - q B (j" 2 + m" 2 ) + it] 
x V,^/w'(q)V> m » ,/,,,< (pi - q - p 3 ) 
+ . . . 

3, 



f d'q 

1 



h(pi) + h(p 2 )o - ep,- q - fp 2+q + 2^ - q B (j" 2 + m" 2 ) + it] 
x [{jm\T = 0}(T = 0\j"m")(j"m"\r = 0}(T = 0|/m') 
x Vo(q)V (pi - q - p 3 ) 

+ (jm\T = OXT = 0\j"m")(j"m"\r = 2}(T = 2|/m'> 
x Vo(q)V 2 (pi - q - p 3 ) 

+ (jm\T = 2){T = 2\j"m")(j"m"\r = 0)(T = 0|/m'> 
x V 2 (q)Vo(pi - q - p 3 ) 

+ (jm\T = 2){T = 2\j"m")(j"m"\T = 2){T = 2\fm') 
x V 2 (q)V 2 (pi - q - p 3 )] 

+ . . . (39) 



Here, the second identity in Eq. (1391 is obtained by using Eq. ( 1261 1 for Gj(p) in the integration of 

G° r (pi - q)Gl t „(p 2 + q) with respect to q Q . 

In spinor BECs, the stable quantum phase of the ground state is determined as the competi- 
tion between the spin-exchange interatomic interaction and the coupling of atoms to an external 
magnetic field via the quadratic Zeeman energy, and thus, the quadratic Zeeman energy usually 
has the same order of magnitude as the spin-exchange interaction: q B ~ \c\\n <K cqh. It can be 
shown that for such an external magnetic field, the spin dependence of intermediate states via the 
quadratic Zeeman energies qsij" +m") in the denominator of the right-hand side of Eq. (l39l only 
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Figure 4: T-matrix of a two-body scattering. Two atoms with momenta Sp3, 7lp4 and magnetic quantum numbers /, m' 
collide to form two atoms with momenta Hp 1 , hp2 and magnetic quantum numbers j, m. The T-matrix is defined as the 
sum of an infinite number of ladder-type diagrams which describe virtual multiple-scattering processes [see Eq. j39H . 



gives a small difference that is negligible up to the order of magnitude we are considering in this 



paper (see Appendix A 1. Consequently, as a good approximation we can take the summation 

2 \8 h ><Sh\ = 1 (40) 

out of the integral. Inside the integral, by using the fact that the T — and T -2 spin channels 
are orthogonal to each other: {T = Q\T — 2) = 0, the T-matrix can be rewritten as 



r jm j m >(pi,p2;p3,p4) ={j,m\T = 0}{T = 0\j'm')To(pi,p 2 ;p3,p4) 
+ (jm\r = 2){T = 2\fm')T 2 (p u p 2 ;pi,p A ), 

where IV(pi, p 2 ; P3, /? 4 ) is the T-matrix in the T spin channel given by 

C d^ 

Tr(puPi;P3,P4) =Vr(Pi - Ps) + ^ I 77^7 G°(j>i - q)G°(p 2 + q)V r (q)V r (pi - q - 
+ 



(41) 

P3) 

(42) 



Here, G°(p) = 1 1 {pa - £p + p + irfi is the spinless non-interacting Green's function. 

The T-matrix Tf(p x ,p 2 ; pi,, /? 4 ) can be expressed in terms of the vacuum scattering amplitude 
for the spin channel T — and 2 as follows (see Appendix A I lfll [35ll : 



Tr(PuP2-,P3,P4) =T r (p,p',P) 

=/r(P,p')+ f ^ 
1 



Vr(p,q) 



hP -^+2p-^ + iT 1 



hrq- 
M 



h 2 a' 2 
M- l1 l 



/r(p'.q). 



(43) 



where -Mfir(p, p') / {4nh 2 ) is the vacuum scattering amplitude of the two-body collision in which 
the relative momentum changes from Hp' to hp. As seen in Eq. (1431 . it can be shown that 
ryQ?i > P2', P3,P4) depends only on the four-vector total momentum hP = hp\ + hp 2 = Hp3 + frp4 
and the relative momenta hp = (hp\ - Hp 2 )/2, hp' = (Hp3 - Hp4)l2, and depends on neither 
Po = [(Pi)o - (P2)o\ /2 nor p' Q = [(p 3 ) - (p 4 )o] 12 (see |Appendix A| >. 
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3. First-order approximation-Bogoliubov theory 



In the approximation to first order in the interatomic interaction, we can neglect the q-integral 
in Eq. d43l because it give a contribution to second order. Indeed, its contribution is smaller in 

magnitude than the first-order contribution by a factor of the diluteness parameter Jnci^r 1, 

where ay- is the s-wave scattering length in spin channel T (= 0, 2) (see Sec. |4). On the other 
hand, in the low-energy regime |p| <K 1 ja^, the momentum dependence of the vacuum scattering 
amplitudes is negligible, and ff(\>, q) reduces to ff = \nfi}afjM in the limit of zero momenta: 
p, q — > 0. The T-matrix then becomes 



(p.p'.P) - (j,m\T = 0)(T = 0|/,m')/ + (j,m\r = 2){T = 2\f,m')f 2 . 



1 jmj'm' 

By using the following relations 

\T = 0){T = 0| + \T = 2)(T = 2| =1, 
-2\r = 0){T = 0| + \T = 2)(T = 2| =F ■ F, 

the T-matrix can be rewritten in the following form: 

Tj lnJm ,(j),p',P) * c 6 jr 6 mm > +c\ ^(F a )jj>(F a )„ 



(44) 



(45) 
(46) 



(47) 



where (F a ) (a = x,y,z) are the components of the spin-1 matrix vector 

1 0^ 



F "T2 



(0 
1 1 
1 0J 



-1 


1 




(1 





) 





-1 


,F Z = 











1 


0; 




lo 





-ij 



(48) 



and Co and c\ are the coefficients of the spin-conserving and spin-exchange interactions, respec- 
tively. They are related to the s-wave scattering lengths as follows: 



c = 



c\ =- 



fo + 2/2 4nH 2 ao + 2a 2 
3 ~ ~M 3 ' 
h ~ /() 4^ft 2 a 2 



M 



(49) 
(50) 



For a convenience, we define a characteristic length scale a (a) that corresponds to the spin- 
conserving interaction in the T-matrix given by Eq. (l4Tt : 



a _ ao + 2a 2 
4^ = 3 ' 



(51) 



from which we have co = Antra /M = h 2 a/M. 

Now, we consider two cases in which the ground state is in the ferromagnetic and the polar 
phase. 

3.1. Ferromagnetic phase 

In the ferromagnetic phase, all condensed particles occupy the j = 1 magnetic sub-level, and 
the condensate's spinor is 



(fi,fo,f-i) = (1,0,0). 
14 



(52) 



The proper self-energies and chemical potential in the first-order approximation, which are illus- 
trated by diagrams in Fig. are then given by 

m^ip) = r A/1 (p/2, p/2, P ) + r wl (p/2, -p/2, P ) 

- Conoid jf +6j t \6y,l) + C\Hq ^ [(•Fa)//>(^'<*)ll + CF w«)lj-] 

= c n (6 jf + 6 M 6 + Cin (J8jj> + 5 M 6 /,i + 5 jfi 5 /, ), (53a) 

= /a*), {p) = r^ vl ( P) o,o) 

C no<5y,l^/,l + Ci« ^(^a);,!^)/,! 

a; 

= c no£/,i£/,i + ci«o£;,i<5/',i, (53b) 

p=r lul (o,o,o) + ^ 

-c h + Cin ^(/v)iiCF ff )n + <? B 

a 

= (c + c l )n a + q B . (53c) 



Here, the quadratic Zeeman energy qg is added to the right-hand side of Eq. (I53cb for the chem- 
ical potential to account for the fact that the condensate is in the magnetic sublevel 7 = 1, whose 
energy is raised by qs due to the quadratic Zeeman effect. The matrix elements of £(p) in Eq. ( l28l l 
are then given by 

KL\\{p)=2{c + c l )n , (54a) 
KE ^(p)=(c + c 1 )no, (54b) 
hl}_\ _ l (p)= (c -cOno, (54c) 

= fiSfl(p) = (c + d)no, (544) 
others = 0. (54e) 

By substituting Eqs ( I53cb and (154-b into Eq. d30b , we obtain the first-order Green's functions: 

po + £/h+(c + c\)n /H 

Gu(p) = ~ j ; , (55a) 

Po ~ W I,p + lt l 

Goo(p) = - (55b) 

PO ~ d>0,p + IT] 

G-i.-i(p) = - (55c) 

PO - W_i jP + IT} 

G n {p) = ^ = _ fo + cOno/^ (55d) 
Po - W l.p + "7 

others = 0. (55e) 
The energy spectra of the elementary excitations, which are given by the poles of the Green's 
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(a) 





(c) 



(d) 
M 



Figure 5: First-order contributions to the proper self-energies (a)-(c) and the chemical potential (d). Here, the squares 
represent the T-matrices while particles belonging to the condensate are not explicitly shown. In fact, in (a), there are one 
condensed atom moving in and one condensed atom moving out; in (b) and (c), there are two condensed atoms moving 
in and two condensed atoms moving out, respectively; in (d), all four atoms are condensed atoms. This convention helps 
simplify the second-order diagrams in Sec. [4] 



functions, are 



fuo Up = -^ep[e° + 2(co + Ci)n ], (56a) 
Hat , p = - qB, (56b) 
fow-i.p - e p ~ 2c\riQ. (56c) 

Thus, the Green's function approach gives the Bogoliubov energy spectra of elementary excita- 
tions as the first-order results Q34J]. 

It will be useful for the second-order calculation in Sec. 0to rewrite the first-order Green's 
functions in Eq. d55l i as follows: 

A.\ n B\ 

Gu(p) = : (57a) 

Po - ^i.p + it] p + a>i P - i7] 

GUip) = Gj\=- C,, p ( —± 1 , (57b) 

\P0- Wl, p + IT] po + Ci>l. p - IT] ) 



where 



ha>i tP + e° + (co + ci)« -fuo hp + e° + (c + ci)«o 

^(co^o 
Lfux>\ v 
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3.2. Polar phase 

In the polar phase, all condensed particles occupy the j — magnetic sublevel, and the 
condensate's spinor is 

(fi, ft, ?-i) = (0,1,0). (60) 
From Fig. [5] the first-order proper self-energies and chemical potential are given by 

HL) l f (p) = r i0 ,/o(p/2, p/2, p) + r , ;/0 (p/2, -p/2, p) 

- c n (Sjj> +Sj t0 Sffl) + ci«o ^ [(^)>./(^»)o,o + (^ a ) j( ;o(^)o,/] 

= c n (5jj +6jfl8ffl) + ci«o(<5/,i<5/,i + <5/-i<S;<-i), (61a) 
RE^(p) = ft2^(p)=r J r / .oo(p,0,0) 

- cono<5;,o<5/,o + ci«o V^Jio^^/.o 

= c no5j,o^/.o + cin (6 jt iSf-i +8j-iSf t i), (61b) 
/x=r 00 ,oo(0,0,0) 

-c «o, (61c) 

The matrix elements of E(p) in Eq. d34l i are then given by 

tSL\\(p) = REl 1 ! _!(p) = (co + ci)«o, (62a) 

t&wip) = 2c n , (62b) 

^^(p) = ^i 2 u (p) = Kl-i(p) = ^ <x (p) = cmo, (62c) 

hX 1 *(p) = Hl^(p)=c no, (62d) 

others = 0. (62e) 



Substituting Eqs. ( 16 let and ( I62l i into Eqs. ( 1361 1. we obtain the first-order Green's functions as 
follows: 

po + (e° + ci« + q B ) Ifr 

Gn(p) = G^(P) = y -f , (63a) 

Po - W I, P + i? 7 

po + (e„ +c «o)/^ 

Goo(p) = r— • ( 63b > 

Po - w o, P + «7 

Gj^Cp) = Gi 2 ,(p) = Gf.^p) = Gf_M = ~ 2 Cin i H . , (63c) 

Po - W I,p + "? 

,~12, \ r,2\, s C n Q /h 

G 00(P) = G Q0(P) = - — 2 ~' (63d) 

Po - W 5, P + "I 
others = 0. (63e) 

The energy spectra of the elementary excitations, which are given by the poles of the Green's 
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functions, are 



hco hp - ?ki»_i jP = J(ejj + q B )(£p + q B + 2ci« ), (64a) 

ftw , P = A /e°(e^ + 2c n ). (64b) 

Here, there is a two-fold degeneracy in the energy spectra: io\ v = w_i, p for the polar phase 
due to the symmetry between two magnetic sublevels j = ±1. Similarly to the ferromagnetic 
phase, by using the Green's function approach, we have obtained the Bogoliubov energy spectra 
of elementary excitations for the polar phase as the first-order results [34]. 

The first-order Green's functions given by Eq. (l63~t can be rewritten in the following form: 



where 



A\ n B\ n 

G 11 (/7)=G_ 1 ,_i(p)= ^ — : ^ — - (65a) 

PO ~ Wl, p + IT] p + L0 Up - IT] 
A 1 n Bq n 

Goo(p) = : (65b) 

Po - w ,p + it] po + w ,p - iri 

G\%(p) =Gi 2 u (p) = Gl_,{p) = G^iCp) 

= - CiJ l - : -1 (65c) 

\fo-wi, p + !77 p Q + a> hp - it] ) 

Goo(P) = G^ip) = -CoJ : l - -1 (65d) 

\p - LL>Q,p + IT] PQ + W ,p - IT] ) 



ho)\, v + 6p + ci« + q B -kioi, p + 6p + ci« + qB 

,P 2ha>ip ' ,p 27wi>i. p 

A ,p = ^ , B Q ,p = P — , (67) 

Cj = ci»o ^ = CQWQ 
' p 2ha>i p ' ' p 2?i<yo,p 

These expressions will be used in the following sections. 



4. Second-order approximation-Beliaev theory 

We now investigate how the effect of quantum depletion at absolute zero alters the energy 
spectra of elementary excitations in an F = 1 spinor condensate of 87 Rb by calculating the energy 
spectra to the second-order in interaction. The spin-exchange interaction for 87 Rb atoms is known 
to be ferromagnetic (ci < 0). Here, we only consider the case of q B < and q B > 2\ci\n for the 
respective ferromagnetic and polar phases, where the corresponding first-order energy spectra 
of elementary excitations show that the system is dynamically stable [see Eqs. d56b and d64l)l. 
On the other hand, when considering the second-order corrections to the first-order results, we 
only need to take into account the spin-conserving interaction since the spin-exchange interaction 
would make a much smaller contribution to the already very small second-order quantities. This 
is due to the large ratio of spin-conserving to spin-exchange interactions of 87 Rb atoms: co/|ci| - 
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200. However, for a usual atomic density in experiments of ultracold atoms, the second-order 
contribution to the proper self-energies from the spin-conserving interaction is of the order of 
con V/ia 3 ~ 0.01co«, which is of the same order of magnitude as the first-order contribution 
from the spin-exchange interaction ~ \c\\n. We may thus expect an interplay between quantum 
depletion and spinor physics. 

4.1. Second-order proper self-energies and chemical potential 

The second-order correction of the proper self-energies and chemical potential involves two 
terms. One is the second-order correction to T<f =02 {p\ , Pi\ Pz,Pa) in the first-order diagrams (see 
Fig- [5j, that is, the q-integrals and the imaginary part of fr=o,i(p, P') in Eq. (143V The other is the 
contribution from the second-order diagrams given in Figs. |6]|9] 



4.1.1. Ferromagnetic phase 

First, we consider the second-order corrections to the self -energies and chemical potential 
that result from the correction to the T-matrix in the first-order diagrams. They are obtained by 
substituting the q-integrals and the imaginary part of frip, p') in Eq. d43l into the first lines of 



Eqs. (I53ab - (l53cb (for more details, see Appendix B 



KL)l(p) : i Im{c (p/2, p/2)}n 6 jr + i Im{c (p/2, -p/2)}n S M S r , 



3 
1 



+ "o 



2\ 



f i!s 

J (2ff) 



7ffij 2 ' 21 Q?) : "o 



ju : n 



o _ c o _ J) 

2\ 



/ f2 



y \ hp + 2[(c Q + ci)n + q B ] - e , - e ^ + ir] 



J (2n) 
J (2n) 



) 3 \ 2[(c + c x )« + q B ] ~ 2el + it] 2$ 



1 



1 



r) 3 \2[(c + ci)n + q B ] - 2e° + it] 2e% >' 
where Im denotes the imaginary part of a complex number, k = q - p, and 

,„„ , . ,~ _ /o(P/2, ±P/2) + 2/ 2 (p/2, ±p/2) 
c (p/2,±p/2) = . 



(69a) 
(69b) 
(69c) 

(70) 



Using the optical theorem for scattering, the imaginary part of an on-shell vacuum scattering 
amplitude fr(p, p') with |p| = |p'| is given by l35ll 



im{/ r (p,p')} = -^ j ^3/r(P,q)/r(p'.q)«y(p 2 -q : 

j dQ q / r (p,q)/^(p',q), 



I6n 2 h 2 



(71) 



where f2 q denotes the solid angle of the on-shell momentum q: |q| = |p| = |p'|. Consequently, 
the imaginary parts of fjr(p/2, +p/2) and co(p/2, +p/2) in Eqs. ( |69l l and d70l i are given in the 
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(al) (a2) (a3) (a4) 




(dl) (d2) (d3) (d4) 




(el) (e2) (fl) (f2) 




Figure 6: Second-order diagrams for the proper self-energy Zjj, (p). The intermediate propagators are classified into 
three different groups, depending on the number of noncondensed atoms moving into and out of the condensate. They 
are represented by curves with one arrow ( — >), two out-arrows (< — >), and two in-arrows (—><—), and are described 
respectively by the first-order normal Green's function Gjy (p) and two anomalous Green's functions GjJ, (p) and G^j, (p), 
respectively. Here, the two horizontal dashes in diagrams (el) and (e2) represent the fact that we need to subtract from 
these diagrams terms containing non-interacting Green's functions to avoid double counting of the contributions that 
have already been taken into account by the definition of the T-matrix and the first-order diagrams [see Eqs. ( IC. 181 and 
IC.19H . As in Fig. [5] here, we use the convention that the condensed particles in diagrams (al)-(e2) are not explicitly 
shown. 
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(al) (a2) (a3) (a4) 




(bl) (b2) (b3) (b4) 




Figure 7: Second-order diagrams for the proper self-energy l}^,(p). Similar to the horizontal dashes in diagrams (el) 
and (e2) of Fig. [6] the vertical dash in diagram (e) represent the fact that we need to subtract from this diagram a term 
containing non-interacting Green's functions to avoid double counting of the contribution that has already been taken 
into account by the definition of the T-matrix and the first-order diagrams [see Eq. jC.40H . 
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(a2) 



(a3) 



(a4) 











r 




y. J 




















.1 





q-p 



(d2) (e) 

p-q f q^-^-q 




P h' 



q *2+->.-q 



Figure 8: Second-order diagrams for the proper self-energy X ,(p). 



(al) 



(a2) 



(b) 



Figure 9: Second-order diagrams for the chemical potential /i. 

22 



second-order approximation as follows: 



Im{/ r (p/2,±p/2)}=-^/ r , 



Im{c (p/2,±p/2)} 



8tt?i 2 



'/ 2 + 2£ 



(72) 
(73) 



where we have replaced /?=-(p, p') on the right-hand side of Eq. ( fTTT i by its zero-momentum limit 

Next, we calculate the second-order contributions to the proper self-energies and chemical 
potential from the second-order diagrams illustrated in Figs.[6]|9]by using the first-order Green's 
functions given in Eq. ( |57| ) (for more details, see Appendix C.l| l. By summing up the second- 
order corrections that arise from the correction to the T-matrix [Eq. d69| il and the contributions 
from the second-order diagrams [Eqs. (IC.2) - (IC.40) 1, we obtain the second-order self-energies 
and chemical potential as follows: 



47Tft 2 



/o 2 +2/ 2 2 



2 \ 



+ 2« 
1 



f- 



q_ 

) 3 



hp + 2[(c + ci)n + q B ] - e q - + ^ e p ~ e q ~ e k + "7/ 

d 3 q /2 {A u , Bi >k ) + 4Ci, q Ci, k - 4 (A u , Ci, k ) + 2A 1 , q A 1 , k 



•n Co 



J (2tt) 3 



h\P0 - Wi, q - Wi, k ) + IJ] 

2 {A,, q , + 4C u Ci, k - 4 (B ljq , Ci, k ) + 2B 1>q B 1;k 



/i(/7 +Wl >q + Wl ik )-/77 

1 2c 



+ 2(c + ci)« + ;/7 



r d 3 ( 

J (2* 



3 fil.q, 



(74) 



^ + 2/^ 



+ "0 

1 



2\ 



r d 3 < 

J (2* 



q_ 

■) 3 



i 



ftpo + 2[(c + ci)n + q B ] - e q - e k + "7 e p ' ~ e q " e k + 
q / Ai k + Bi jk - 2Ci, k 



2 f d 3 q 

" 0C() J W 



) 3 \?i(po - w , q - <wi,k) + ;'?7 



- e q - e k + 2(c + ci)« + <7b + »7 



+ c 



r_d^ 

J (2^ 



(75) 
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-i\p\Mn 



+ «o 



'JS + 2JZ 



r d 3 ( 

J On 



) 3 



7yj + 2[(c + Ci)« + 9b] - fq ~~ e k + "7 e P ~ 6 q ~ e k + "7 



f _d*q 
J (2ff) 



A, k + B 



l.k ■ 



2C 



l.k 



) 3 \h(po - - wi,k) + «7 



1 



Co 



f B 



(76) 



< m (p)=«E^(p) 

f/o+2/l 



f_d^q 

J (2jt) 



) 3 \2[(c + a)n + q B ] ~ 2e° + b] 2e\ , 
+ m>c§ J ^ (2{A 1 , q ,B 1 , k }+6Ci, q C lik -2{A 1 , q + B 1 , q ,Ci,k}) 
1 1 

f d 3 ( 

J (2tt: 



+ c 



q_ 

T) 3 



Ci,,+ 



h[p + u> Uq + wi, k ) - it] 

CQHQ 



, 3 JJ (2»r 



2c 



C d 3 q f d 3 ( 



2e° - 2(c + ci)« - 07/ 

1 1 
2[(c + c l )n + q B ]-2e® + iri 2e° 

c «o 



-Ci, q + 



2e q - 2(c + ci)n - »7 



(77) 



(78) 



where k = q - p and [A h(l , Bi >k } = Ai iq B ljk + Ai, k Bi >q . 

Here we consideronly the case in which the external magnetic field satisfies q B ~ \c\\n «: con 
(see Sec. 12.31 ). and ignore terms smaller than co« Vna 3 , which is the order of magnitude of the 
second-order approximation under consideration. Then Eqs. d74b . d77b . and d78l l for '(/?), 



5ii^(p)» ar, d J u ^ ) ' respectively, are the same as those for a spinless Bose-Einstein condensate 



,(2) 



lj. It is because the condensate is in the j = I sublevel, and the elementary excitation given 
by Sj[' 12(2) (/?) is the density-wave excitation as in a spinless system. Consequently, it has a 
phonon-like second-order energy spectrum in the low-momentum regime (eJJ «: cqu): 



/ 7 _ v r 3 _ (e«) 5/2 

Tlpo = y + Vnofl 3 j V 2 "o(co + ci)yep - '^jjj^oco V"o« 3 - ^ 5/2 

5/2 



1 + V"oa 3 ) V2«o(co + ci) Je£ - i^^-n c y/n^a 
3 V 71 ' / y 8U 

(0 



(4) 



1 + 



j= Vna 3 j y/2n(c + ci) ^ - i 



3 
~80 



(n co) 5/2 

5/2 



(nc ) 



5/2' 



(79) 
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where a and a are defined in Eq. ( BTT i. Here, in the derivation of the last line of Eq. d79l l. we used 
the expression for the condensate fraction in a homogeneous system 



«o = n 1 — 



3 



(80) 



The first term (real part) on the right-hand side of Eq. ((79} shows an increase in the sound velocity 
of a density-wave excitation due to quantum depletion, while the second term (imaginary part) 
is the so-called Beliaev damping, which shows a finite lifetime of phonons due to their collisions 
with the condensate. The second-order contribution to the chemical potential is given by |l|] 



fj. (2) = "oco V«o5 3 - 



(81) 



Now, to evaluate XLi^(/?) we take a Taylor expansion of it around po = ol>o p , where Ha>o tP is 
the first-order energy spectrum given by Eq. (156b) : 



OP) 



yll(2)/ n 

2 oo (Po 



W0,p) + 



d yU(2) 
OZj 00 



(p) 



dpo 



(po - w .p) +0 [p - w ,p) 



(82) 



We can stop at the linear term in this Taylor expansion, provided that the difference between the 
second-order energy spectrum and the first-order one is small: 



\po - W ,p| <sc 



v ll(2), 
V \P0 



W ,p) 



CQ»Q 

h ' 



(83) 



[dl^/dpo] (po = w .p) 
which is justified by the fact that the system is a dilute weakly interacting Bose gas, and will be 



confirmed later by the second-order energy spectrum obtained below in Eq. (II 19l l. This will be 
discussed in more detail at the end of Sec. 14.2.11 

It can be shown that the imaginary parts of E^ 2 \po = u>o,p) and [dE^' 2 ' /<9po] (Po = k>o,p) 
vanish for any value of p (see Appendix D. 1 1, which results in 



Irr£ fV) = + tf[(A)-^o, P ) 2 l 



(84) 



This result implies that there is no damping for the elementary excitation given by up to 

the order of magnitude under consideration. 

For the real parts of Z^^ipo = a>o, P ) and [<9£ ( ' <2, /dpo] (po = coo,p), we can make their Taylor 
expansions around p = in the low-momentum regime el] <sc cq«q: 



Re^ (2) (/7 



dR^ i2 \p) 



ajo.p) =ReI. l Q Q (2 \p 



dRel,l^ 2) (po = w ,p) 



<^0,p) 



2 



l P =o 

W0,p) 



d(cj up ) 2 



p=0 



do 



x «T, P + 



p=0 



dpo 



dReZ l f\p) 



dpo 



x E Ux> + ~ 



d ldR^lf\p) 



Po=«o. P ,P=0 

d 2 



dco hp \ dp 



PO-^O.p 
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X OJ 



p=0 



Up 



P =o 



(85a) 



(85b) 



where wi p is given in Eq. d56al >. Note that p = is equivalent to wi jP = 0, and e° <K co«o is 
equivalent to hu>\# <sc co«o- 

With straightforward calculations, we obtain (see |Appendix E.r| for details): 
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{ dp 


1 <9 2 


fdR^ 2 \p) 




{ dp 
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Po ~ w 0,p 
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p=0 
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V«o5 3 



(noco) 2 



(86) 
(87) 

(88) 

(89) 

(90) 
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Hence, we obtain: 
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1 - 



49 / Swi, p 



+ [q? - w , p ) 2 ] . 



1200 \ noco 
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40 
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"oco 
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(92) 



Equation ( f92b shows the modification of the self-energy £ii(p) due to the effect of quantum 
depletion. The first term in the first line is the value for p = 0, po = 0, the second term in the 
first line is the correction for a nonzero momentum, while the second line is the correction for a 
nonzero energy. It can be seen that the self-energy E^fp), which describe the effect of interaction 
with other particles on the propagation of a quasiparticle, decreases with increasing momentum 
or frequency. 



4.1.2. Polar phase 

Following a procedure similar to the ferromagnetic case, the second-order corrections to the 
self-energies and chemical potential that result from the correction to the T-matrix in the first- 
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order diagrams are given by 



KL)\(p) :/Im{c (p/2,p/2)}«o^ 7 + *Im{e (p/2, -p/2)}n ^, S/,o 



no 



fS + 2 fi 

3 
1 
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J (2* 



hZf/Hp) : no 
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p fc q fc k 



) 3 \ hp + 2c n - eq - + ('77 



J (2»r) 
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) 3 \2c «o-2e5 + /77 2e° 



5 fyofy.o, 



) 3 \2c «o-2e^ + /77 2e£ 



(93a) 
(93b) 
(93c) 



By summing up the second-order corrections that arise from the correction to the T-matrix 
[Eq. j93l l and the contributions from the second-order diagrams (see |Appendix C.2\ , we ob- 
tain the second-order self-energies and chemical potential as follows: 
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(94) 
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1 
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ft(p - ^i.q - 0»l,k) + it] 



Co, q + 



C()«0 
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(95) 



(96) 
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J (2?r 



2e q - 2c n Q - irj 

I (2^ ( 2c o«„-2 e o + /, + 2^) + 2C ° / (2^ ^ + 

r \ 



(97) 



_ Co, q - 



CQ«Q 

2c «o - 2eq + irj 
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(98) 



It can be seen that the integrand of the q-integral in each of Eqs. d95l l, (l97l i. and d98l l is a sum of a 
term that contain only Ao, q; k, fio, q ;k, Q),q;k, wo,q ; k, co«o an d a term that contain only Ai q; k, fii, q; k, 
Ci, q; k, Wi, q; k, Ci«. By rewritting the corresponding q-integrals using dimensionless variables 
eq/(co«o) and e q /(ki|no), we find that the value of the latter integral is smaller than that of the 
former one by a factor of VkiT/co 1, and thus, the latter integral can be ignored. Here, 
we used hpo fiwo.p <K Vki|«oco«o <K Wo for the low-momentum region ep* <K |ci|« under 
consideration for the case of the polar phase. Consequently, (p) and // 2) are the same as 

the second-order self-energies and chemical potential of a spinless Bose-Einstein condensate U]. 
Namely, the second-order contribution to the chemical potential is given by 



P- 



(2) 



3tt 2 



n c 



\/«o5 3 - 



(99) 



1 1 ' 17.(2) 

Here, the elementary excitation given by S Q0 ' (/?) is a density-wave excitation as in a spinless 
system. It, therefore, has a phonon-like second-order energy spectrum in the low-momentum 
regime: 



5/2 



(nc ) 



5/2' 



(100) 



On the other hand, in Eq. (|96] l for Xj 2 *; 2 ', the factor ci«o, which arises from Ci iq , can be taken 

out of the integral, and thus, Ej _j is negligibly small compared to the order of magnitude under 
consideration: 



= + 



l"o J«o^ 



(101) 



Now, to calculate £[j (2) (/?) we make its Taylor expansion around po = where Swi iP is 
the first-order energy spectrum given by Eq. (164a) : 



-11(2)/ 



(p) = m a) (po = ^i,p) + 



dpo 



(po - ui,p) + (po - wi.p) 



po=<yi,p 



(102) 



We can stop at the linear term in this Taylor expansion, provided that the difference between the 
second-order energy spectrum and the first-order one is small: 



v ll(2), \ 
gjj AgO - ^l,p) 

[3Sj| (2) /3p ](po = wi, p ) 



Po - wi, p «: 



(103) 



which is justified by the fact that the system is a dilute weakly interacting Bose gas, and will be 
confirmed later by the second-order energy spectrum obtained below in Eq. (1140b . 

It can be shown that the imaginary parts of 2[[ <2) (po = wi,p) and [<9X[[ (2) /<9f>o] (po = <±>i,p) 
vanish for any value of p (see|Appendix D.2|), resulting in 



imz;; (2) (p) = o+o[(po-wi. P ) 2 ] 
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(104) 



Physically, this implies that there is no damping for this elementary excitation up to the order of 
magnitude under consideration. 

For the real parts of S[J (2) (po = wi iP ) and [3£|[ (2) /(3/?o] (po = <^\,p), we can make their Taylor 
expansions around p = in the low-momentum regime <K \c\\no cono: 
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(105b) 



where wo, P is given in Eq. d64bb . Note that p = is equivalent to «o,p = 0, and e° <K co«o is 
equivalent to fiwo.p <sc co«o- 



With straightforward calculations, we obtain (see Appendix E.2 for details): 
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Therefore, we have 



HReI.\\ i2) (p) = -j-2 n o c o V«o5 3 
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(112) 



Equation (II 12b shows the modification of the self-energy due to the effect of quantum 

depletion. The first term in the first line is the value for p = 0, po = 0, the second term in 
the first line is the correction for a nonzero momentum, while the second line is the correction 
for a nonzero frequency. It can be seen that because (q B + Ci«o)/ V?b(9b + 2ci«o) 1 for 
any q B ^ 2|ci|«o> the self-energy Sj [(/?), which describe the effect of interaction with other 
particles on the propagation of a quasiparticle, decreases for increasing momentum or frequency, 
regardless of the strength of the external magnetic field. Similarly, we have 
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(113) 



4.2. Second-order energy spectra of elementary excitations 

With the second-order self energies and chemical potential obtained in Sec. 14. II we are now 
in a position to evaluate the second-order energy spectra of elementary excitations, which can 
be obtained from the poles of the second-order Green's functions. As shown in Sec. 14.11 there 

1 1*12 11*12 

is always one density-wave elementary excitation, which is given by G\\-^ and for the 

ferromagnetic and polar phases, respectively. It has a linear dispersion relation as the phonon 
mode in spinless BECs [see Eqs. d79l ) and ( 11001 )1. As a consequence of quantum depletion, the 
sound velocity increases by a universal factor of 1 + (8/ \Jn) Vna 3 , while there appears the so- 
called Beliaev damping due to the collisons between quasiparticles in the elementary excitation 
and the condensate. The second-order energy spectra of the other elementary excitations will be 
discussed in the following. 
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4.2.7. Ferromagnetic phase 

The poles of the Green's functions Goo(p) and G-\-\(p) given by Eqs. (f30b and OTT ) are the 
solutions of the following equations: 

GooO) : hpo =£p + ^oo(p) 

=4- qB + [^\p)-^} 

=hcj , v + [n^\p)-^], (114a) 
G_i -\{p) : hp =t%-p. + q B + Tffi" _ x (p) 

=haj- Uv + [asif_\(p) - // 2) ] . (1 14b) 

Here, we used the first-order self-energies and chemical potential, which are given in Eq. ( 154) . 
and the first-order energy spectra ha>o--i, P given in Eq. (T56b . Note that on the right-hand sides of 
Eqs. (11 14ab and (11 14bb . the self energies are functions of both po and p. 

By substituting Eqs. ( f84b and d92t for sjj (2) (/j) and Eq. d8TT > for the chemical potential // (2) 
into Eq. (11 14ab . the equation for the pole of GooO?) becomes 

po = o} , p + a p +yS p (po - ^o.p) , (115) 

where o- p , yS p are the lowest-order coefficients in the Taylor expansion around po = wo,p and 
given by 

fia p = - . n c V"()5 3 ( — | , (H6) 



p " 720/r 2 



n c 



*-3?^-i5?^(^) 2 - 



Using the fact that fi p ~ y/tioa 3 «: 1, the solution to Eq. dl 151 ) is given by 



or, 

Po = <^0,p + 



<^0,p + ff P 



1 49 rrj/^i-Pl / 11H n 
iV"oa j • (H8) 



n/zOn 1 \ hqcq 

In the low-momentum region e p ' <K co«o> ^i,p given by Eq. d56ab can be approximated by 

hiQ\ p ^2(cq + ci)no€p. Substituting this and Eq. (I56bb into Eq. ( 1118) . and neglecting all the 
terms that are smaller than the order of magnitude under consideration, we obtain the energy 
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spectrum: 

49 V"()5 3 
720tt 2 n Q c " ^ u ' h '"p 



, o ^ V"oa J o 
%o =£p -qB- „ on 7 x 2c nae„ 



49 



45 



V"oa 3 )ep -t 



*[l--^-^\4-q B . (119) 



45 

Here, we used Eqs. ( fBTt and (fSOb in deriving the last two equalities. Equation ( II 19b shows that 
the energy spectrum of the elementary excitation with spin state j = has a quadratic dispersion 
relation at low momenta, which can be expressed using an effective mass M e ff as 

hpo = ^t s - qB ' (120) 



M 

M eff = — . (121) 

45 -0F 



where 



Compared with the first-order energy spectrum ftwo.p [see Eq. (!56bH . the energy gap remains 
unchanged while the effective mass M e ff of the corresponding quasiparticles increases by a factor 
of 1 /[l —49/(45 y/jf) Vna 3 ]. From Eq. ( l80t . it can be seen that enhancement factor of the effective 
mass is proportional to the number of quantum depleted atoms, both of which are proportional 
to Vna 3 . This can be understood as the effect of the interaction between a quasiparticle and the 
quantum depleted atoms, which hinders the motion of the quasiparticle. 

Furthermore, because the imaginary part of ^\^ T> vanishes up to the order of «oCo y/noa 3 [see 
Eq. d84lil. the damping of this spin-wave elementary excitation is much smaller than that of the 
density-wave excitation mode (the mode with spin state j = 1). In other words, the lifetime of 
the corresponding magnons is much longer than that of phonons. (The fact that the excitation 
with spin state j = and its quasiparticles can be identified as a spin wave and magnons will be 
discussed in Sec.|5]below.) This agrees with the mechanism of the Beliaev damping via collisions 
between quasiparticles and condensed atoms. Physically, the Beliaev damping can be understood 
by considering the conservation of momentum and energy in the collisional process between a 
magnon and a condensed atom. Because Co/|ci] - 200 » 1, the interaction between 87 Rb 
atoms in a scattering process is dominated by the spin-conserving interaction. Consequently, the 
collision between a magnon (spin state j = 0) and a condensed atom (spin state j = 1) would 
produce another magnon (spin state j = 0) and a phonon (spin state j = 1). This is illustrated in 
Fig. [TO] 

The conservation of the total momentum and energy in the collision requires that the follow- 
ing simultaneous equations be satisfied: 

Pl=P3+P4, (122) 

ft 2 p 2 ?l 2 p 2 

^ + 2M e t = -^ + 2M e ; + ?lVslP41 - (123) 
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magnon p 3 ^^r magnon 



P2=°.-T -^P 4 

condensed 

i =i j =1 ^ pnonon 

atom J 2 1 >* K 

Figure 10: A collision between a magnon in spin state j = and a condensed atom in spin state j = 1. Because the 
dominant interaction is spin-conserving, the collision produces another magnon in spin state 7 = and a phonon in spin 
state 7=1. The condensed atom is represented by a dashed line. 



Here, we used the following energy spectra of magnons and phonons: 



+ (124) 

r pho 



E^=hvM (125) 
where the effective mass M e g and the sound velocity v s are given by [see Eqs. dl 19l l and (T79i >l 

M eff = , M , , (126) 
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45 Vj? 



, . (go + cQn 

1h — — yna? A / . (127) 



From Eqs. dl22b and (1123b . we obtain an equation 

fi-IPl +P3|C0Sf? 



2M t 



eff 



= v s , (128) 



where 8 is the angle between pi + P3 and pi - p3 = p 4 . By using Eqs. ( 11261 ) and ( 11271 ) for the 
effective mass and sound velocity, we can evaluate the ratio of the left-hand side to the right-hand 
side ofEq. (fl28l : 

^IPl +P3||COS0| ^ H\pi\ 



2M eff V s M eff v s 



\/2Mc « 



M eff 



1_ _J8_ V^3 



45V? 2C: 



-0 



1 + 4= V«fl 3 V C + C! 

-0(1). (129) 
Here, we used | j>3 1 < |pi|, which results from Eq. (1123b . and |pi| <K yJlMc^ii for the low- 



momentum region e ' <K cqii where the obtained second-order energy spectrum [Eq. (11 19H is 
valid. From Eq. ( 11291 ), it is obvious that there is no collision between a magnon and a condensed 
atom in which the momentum and energy conservations [Eqs. (11221 ) and ( 11231 ). or equivalently, 
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Eqs. (1 122b and ( 11281 )1 are satisfied. Consequently, the damping of magnons vanishes up to the 
order of magnitude under consideration in this paper as shown by Eq. d84l ). 

Similarly, the energy spectrum of the elementary excitation given by G-\-\(p) at low mo- 
menta e l p <K co«o is given by 

hpo 7kj-i,p - 7 n c V"oa d 

720tt- \ « c / 

n _ 49 -Jnod 3 

e _ - 2ci« - „ OA j x 2(c + ci)« e 

1 HOn 1 hqCq y 



V«oa 3 j e p ' - 



1 " V 



49 \ 

1 -V«fl3 e£-2cm (130) 

45 / P 

It can be seen from Eq. (II 30b that the energy spectrum of the excitation with spin state j = — 1 
also has a quadratic dispersion relation at low momenta, and compared with the first-order energy 
spectrum, the energy gap remains unchanged while the effective mass increases by the same 
factor of 1/[1 - 49/(45 y/n) V«a 3 ] as the excitation with spin state j = 0. 

Now we are in a position to evaluate the validity of the a priori assumption that the difference 
between the second-order and first-order energy spectra is small [see Eq. (l83ll. This assumption 
has been used in Sec. 14. li as the self-energies were Taylor expanded around po = cjq p [Eq. (l82l)l, 
and the expansions were stopped at the linear term. The condition for the validity of the Taylor 
expansion can be obtained from Eq. (II 15b , that is, 

/? p |/?o - w ,p| •« <*p. (131) 

By using Eq. ( II 181 ) for the second-order energy spectrum, we find that the left-hand side of 
Eq. ( 11311 ) is almost equal to a p j3 p , and thus, Eq. ( 11311 ) is satisfied, provided that 



j6 p ~ « 1. (132) 

Equation (1132b is nothing but the diluteness condition, and it is usually satisfied in conventional 
experiments of ultracold Bose gases. 

4.2.2. Polar phase 

For the polar phase, there is a two-fold degeneracy in the energy spectra of elementary ex- 
citations due to the symmetry between the j — +1 sublevels. The poles of the Green's function 
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Gu(p) ( or equivalently, G_i -\(p)) given in Eqs. ( f36b and d37b are 
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(133) 



(134) 



2ftwip 

Here, we used the first-order self-energies and chemical potential given in Eq 
energy spectrum given in Eq. (164aK and the fact that Ej 2 ^ and + E_ - 2// 2) /ft are 
smaller than \c\\no/h ~ q B /h ~ oji, p in the low-momentum region e° <K |ci|«o [see Eqs. 
(USD, 053 . <fri2l . and ( fTT3l l. 

By substituting Eqs. (|99]l, dTOTb. ( fT04b . ( fTT2b , and ( fTT3l into Eqs. (fT33l and ( fT34l l, we obtain 
the equation that determines the pole of G\\(p): 
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Using the fact that y8 p ~ -\/n^a 3 <s 1, the solution to Eq. d 1 35b is given by 
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In the low-momentum region ej; <K |ci|«o ~ 1b co«o, a>i p and coo, P can be approximated as 



Hcj Up = J(e° + q B )(£p + q B + 2cj«o) - ^Jq B (q B + 2cm) + (iB + cm) 

V<7b(<7b + 2ci« ) 

Sa» ,p = -\/ e p( e p +2c n ) - -y/2c « 



(139a) 



(139b) 



Substituting Eq. ( 1 1 39b into Eq. ( 1138b . we obtain the energy spectrum which is correct up to the 
second order: 
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yJq B (q B + 2ci«) p 



Here, in deriving the last equality, we used «a 3 «: 1 and 
(?B + cin ) (q B + C]n) 



y/q B (q B + 2cm) y/q B (q B + 2c in) 
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From Eq. (1140b . it can be seen that the energy spectrum of the elementary excitation given by 
G\i(p) has a quadratic dispersion relation at low momenta, which can be expressed using the 
effective mass M e g as 



^2p2 

hpo = ttt — + ylqBiqB + 2ci«). 

ZM eff 

The effective mass depends on the quadratic Zeeman energy q B as: 

y/q B (q B + 2c in) M 



M eff = 



(q B + cin) i 
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(142) 



(143) 



Compared with the first-order energy spectrum given by Eq. (11 39ab , the energy gap remains 
unchanged while the effective mass increases by a factor of 1 / [l - 49/(45 yjn) Vna 3 ] as a conse- 
quence of quantum depletion. It can be seen that the effect of quantum depletion on the effective 
mass is characterized by the same enhancement factor, regardless of whether the system is ferro- 
magnetic or polar, and independent of external parameters of the system. Furthermore, since the 
imaginary part of S[[ (2) vanishes up to the order of «oco s/noa 3 , the damping of this spin-wave 
excitation (see Sec. [5]) is much smaller than that of the density-wave excitation (with spin state 
j = 0). In other words, the lifetime of the corresponding magnons is much longer than that of 
phonons. This agrees with the mechanism of Beliaev damping and can be understood by con- 
sidering the conservation of momentum and energy in a collision between a quasiparticle and a 
condensed atom (see Sec. 14.2. IT ). 
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5. Transverse magnetization and spin wave 



As shown in Secs.[3]and[4] the elementary excitations of a spinor BEC include both density- 
wave and spin-wave excitations with quasiparticles being phonons and magnons, respectively. 
The energy spectrum of phonons with a linear dispersion relation can be experimentally mea- 
sured by using the neutron scattering or the Bragg spectro scop y. The former has been widely 
used in experiments of the superfiuid helium 091 4C , l4l[|42l 14311 . while the later has been applied 



to measurements of ultracold atoms 1144 145 



4611 . Similarly, the neutron scattering has been 



used to measure the dispersion relation of magnons in solid crystals 0471, 14811 . though its applica- 
tion to ultracold atomic systems is limited by a huge difference in energy scales between neutrons 
and atoms. The Bragg spectroscopy can also be generalized to measure the energy spectrum of 
magnons by using appropriately polarized laser beams to make spin-selective transitions. In this 
section, we propose an alternative experimental setup to indirectly get information of the effective 
mass of magnons in a spinor BEC of ultracold atoms. In ultracold atom experiments, atoms can 
be optically excited to a higher energy level to produce an appreciable number of quasiparticles. 
Furthermore, an in-situ, non-destructive measurement of magnetization can be made with a high 
resolution l49ll . It has been applied in a wide range of spinor BECs' experiments to investigate, 
for instance, the formation of spin textures and topological excitations as well as their dynamics 

1 13, m mm Ei- 

First, we show that in a spinor BEC, the elementary excitations with quadratic dispersion 
relations (see Secs.[3]and|4|i can be interpreted as waves of transverse magnetization. The trans- 
verse magnetization density operators F^(r, f) and P y (r, f) in the Heisenberg representation are 
defined by 

P + (r, t) ee FJr, t) + iF y (r, t) = V2 [#(r, #o(r, t) + i%(r, r)#_i(r, t)] , (144a) 
FJr, t) ee FJr, t) - iF y (r, t) = V5 [#J(r, #i(r, t) + ^(r, f)# (r, f)] ■ (144b) 

The squared modulus of the transverse magnetization is written in terms of P+(r, f) and F_(r, t) 
as 

P 2 ± (r,t)=P z x (r,t) + P*(r,t) 

= X - [P + (r, f)F_(r, t) + F_(r, f)^+(r, 0] . (145) 

5.1. Ferromagnetic phase 

The condensate wavefunction is given by 

<p = Vno(l,0,0). (146) 

The lowest-order transverse magnetization P+(r, f) and F_(r, t) are then given by 

F + (r, t) = V2^Ao(r, t) = ^^^^(0, (147a) 

k 

P-(r, t) = V2^%(r, t) = >/2*>5] e "*' ra 0kW' (147b) 

k 

where 2o,k is the Fourier transform of the field operator ^o( r ) = oo(r) as defined in Eq. ( fT5l ). 
Here, we replaced operators (^i(r, t) and ij/\{r, t) in Eqs. ( 11441 ) by the condensate wavefunction 
01 = V«o- 
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At the lowest order (first-order) approximation, by using the Bogoliubov transformation 

b\,k = u k a hk - v k a\ _ k , (148) 
the Hamiltonian for the noncondensate part given by Eq. (JT3J can be effectively written as 

), (149) 

k^O k 

where fi.w ± i ; o,k are the lowest-order energy spectra of elementary excitations given by Eq. ( |56"1 ). 
The system's ground state is defined as the vacuum of annihilation operators b\ k , ciQ k and S-i.k: 

£i*lg> = 0, (150) 

ao, k |g) = o, d5i) 

a-i*lg>=0. (152) 

If a panicle with momentum hko and spin j = is created above the ground state, the system 
is in an excited state given by 



|e> = 4jg>. (153) 



We now calculate the equal-time spatial correlation of transverse magnetization in the system for 
this plane-wave excited state. Using Eq. ( 1 147b . we have 



(e\F + (r, t)P-(r', f)|e> =2« J] e'' (kr ~ k ' r '\e|a o ,k(0< k ,(0|e> 

k,k' 

=2« J] e^-V^e-^*-^' (<5 k , k , + 6 k , ko 6 k>M ) 

k,k' 

=2« J] + 2«oe , ' ko (r " r ' ) , (154) 



k 

where 

ao,k(0 = e^'a , k e-^' = e^ao,*, (155) 

aj k (f) = e^'al k e-^' = e^'a^. (156) 

Here, the first term in the last line of Eq. (1154) , which is proportional to S(r - r'), describes the 
self-correlation, and it exists for all states including the ground state. Similarly, we have 

<eF_(r, t)F + (r', r)|e> =2« £ e"'' (kr " k ' r ' ) (e|aj j£ (0a ,k'(?)|e> 

k,k' 

=2n e"*°' (r " r ' ) . (157) 

Using Eqs. (11541 i and ( 1157b . we obtain the spatial correlation of transverse magnetization in the 
system: 



<e|F + (r, t)P-(r', t) + P_(r, i)P + (r', /)[e> - <g|F + (r, f)F-(r', f) + P_(r, f)^ + (r', 0lg> 



<e|F x (r, • F x (r', f)|e> - <g|F x (r, f) ■ F ± (r', f)lg> 
1 

~2 

=2n cos [ko ■ (r - r')] . (158) 
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Here, we subtract the correlation in the ground state from that in the excited state to remove the 
self-correlation term in Eq. (1154) which is not a physical observable. Equation (1158b shows that 
the elementary excitation given by Eq. (1153l l can be interpreted as a spatial modulation of the 
transverse magnetization, or a transverse spin wave. 

Now let us assume that one particle is excited to create a superposition state of different 
momenta: 



I e™) 



= J d 3 k/(k)< k lg>, 



(159) 



where /(k) is the weight of the superposition. In a manner similar to the above calculation for 
the plane-wave excited state, the expectation value of the squared modulus of the transverse 
magnetization density F^(r, t) with respect to this excited state is given by 



<e sp |Fi(r,f)|e sp )-<g|Fi(r,f)lg> 



1 



<e sp |F + (r, t)F-(r, t) + F_(r, t)F + (r, f)|e sp ) 



- {g\F + (r, t)F4r, t) + F_(r, t)F + (r, f)lg) 



Z 



k.k' 



e /'(k-k')-r e -!(w ,k-Wo.k' )' f ft' 



/(k')7(k) + c.c. 



-2rii 



z 



!(kr-(j , k /) 



/(k) 



(160) 



The expression inside the vertical bars in the last line of Eq. ( 11601 ) is nothing but the time evo- 
lution of a wave packet, which is initially constructed by a superposition of plane waves with 
a weight function /(k). Although the results in this section are derived at the lowest-order ap- 
proximation, as we move to the second-order approximation, the physical properties of the ele- 
mentary excitations do not change (see Sec. [3] and|4]i. Namely, the elementary excitation given 
by Eq. d 153b always has a quadratic dispersion relation and can be interpreted as a transverse 
spin wave. The only difference between the first-order and second-order results is the enhance- 
ment factor for the effective mass of magnons as a consequence of quantum depletion [see, for 
example, Eqs. ( I56bb and dl 19H . Therefore, we can apply the time evolution of the transverse 
magnetization density for a spinor wave packet, which is given by Eq. (11601 ), to the second-order 
approximation with just a replacement of the first-order energy spectrum fi^o.k = Q - 1b by the 

second-order one fkof^ = [l — 49 V«a 3 /(45 V^)] e£ ~ <1b- 

As an example, let us consider a Gaussian wave packet in one dimension, which is a super- 
position state with spectral weight f(k) in momentum space given by 



/(*) = e 



(161) 



This Gaussian wave packet has a width of the order of d in the coordinate space and the center of 
mass moves with momentum Mo. Although a generalization to three dimensions is straightfor- 
ward, a quasi-one-dimensional atomic system is relevant to the experiments of ultracold atoms 
which are tightly confined in the radial direction. Now, we can see how a quadratic dispersion re- 
lation fud^l = h 2 k 2 /(2M e jf) with M e g given by Eq. ( 1121b affects the propagation of a spinor wave 
packet. Note that the energy gap -qe in the energy spectrum has no influence on the squared 
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modulus of the transverse magnetization density because its contribution drops out upon taking 
the absolute value in Eq. j 160b . We then have 



^^m-iS dfcexp 

* -00 

V2^ 



z kx - 



t\-d 2 (k-k y 



v_ 

'27T 



V M2 + it 



exp 



hk 2 

2M eff 

-M eff x 2 - 2id 2 k a (hk t - 2M eff i) 



4<i 2 M eff + 2iHt 



(162) 



and thus, 



2tt 



< 2 V 



2d 2 1 + 



4rf 4 Af 2 , T 



: exp 



-(x - hk Q t/M eff ) 2 



2«f 2 1 + 



(163) 



From Eq. (1163b . it can be seen that due to the nonlinear dispersion relation, the transverse mag- 
netization of a spinor wave packet expands in space during its propagation with a group velocity 
v g = hko/M e ff. The time dependence of the width of the wave packet is also governed by the 
effective mass M e g: 



d(t)=dJl + 



h 2 t 2 
4dMf s 



(164) 



Consequently, by measuring either the group velocity or the rate of expansion of the transverse 
magnetization of a spinor wave packet, we can find the effective mass of magnons. 

To produce a spinor wave packet, which is a localized excited state as given in Eqs. ( 11591 ) and 
(11611 ). a small region of the atomic condensate can be exposed locally to a pair of laser beams 
which couple the states in the ground-state manifold to the electronically excited states. Via a 
Raman optical transition, which is a two-photon process, a fraction of the atoms in the j = 1 
sublevel are transferred locally into the j = state (see Fig. ITTt. As can be seen from Eqs. (l85"b 
and j92l , which are Taylor expansions in powers of the momentum, the second-order energy 
spectra obtained in Sec. |4]are valid only in the low-momentum region ej,' «: co«o- Using the 

■ 87t 



(165) 



parameters of Rb 15511 . the maximum momentum hp m . dx is given by 



V8^ 



10 7 m -1 . 



Therefore, as a superposition of plane waves with momenta in the above region, the spinor wave 
packet should have a width of the order of 



Ax 



1 



» 10~ 7 m. 



(166) 



The condition (1166b turns out to be well satisfied with the parameters of laser beams used in the 
experimental setup. The pair of laser beams is set to be perpendicular to the single laser beam 
used as the trapping potential (see Fig.fT2l. The wavelength of the pair of laser beams that couple 
the ground-state manifold to electronically excited states is of the order of 0.5 fim and their beam 
waist is a couple of the wavelength, i.e., of the order of a micrometer. Finally, using Eq. (1164b 
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Figure 11: Two-photon Raman optical transition used to transfer Rb atoms between the 7=1 and j = spin states 
in the ground-state manifold 5 2 Si/2, F = 1. To produce a localized wave packet of transverse magnetization, atoms in 
a small region of the atomic cloud need to be transferred from the j = 1 to j = spin state for the ferromagnetic phase 
(a), and from the j = to j = 1 spin state for the polar phase (b) (see Sec. 15.21 . Here, <x + (<x~) denotes the right (left) 
circularly polarized light, and n the linearly polarized light. 
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Figure 12: A schematic configuration of laser beams used to produce a wave packet of transverse magnetization. A pair 
of laser beams both along the x-axis are used to transfer atoms between the j = I and 7 = spin states: one is linearly 
polarized (rr) in a direction parallel to the external magnetic field, which determines the quantization axis (z-axis), while 
the other is linearly polarized in a perpendicular direction (v-axis), which can be regarded as a superposition of two 
circular polarizations (<r + + cr~). An additional single laser beam along the z-axis is used as a trapping potential. 
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and the parameters of 87 Rb 15511 . we can estimate the time it takes for the spinor wave packet 
to expand to the entire atomic cloud. For an atomic cloud whose axial length is 100 fim, the 
evolution time of the spinor wave packet is about 40 ms. It is well within the lifetime of the 
condensate, which is of the order of a second. Furthermore, the enhancement of the effective 
mass of magnons as a consequence of quantum depletion is manifested by a difference in the 
width of the spinor wave packet, which is of the order of 1 /mi after 40 ms of propagation. 

5.2. Polar phase 

The condensate wavefunction for the polar phase is 

0= V^o(0, 1,0). (167) 
The lowest-order transverse magnetization F + (r, f) and F_(r, /) are then given by 
F + (r, f) = [#(r, t) + #_,(r, t)] 

k 

= «""** K k (0 + a_i,-k«] , d68a) 

k 

F_(r, f) = V2^ [#!(r, t) + 0-:,(r, t)] 

= ^Yj [e' kr «i,k(0 + e- ikr a[ l k (t)] 

k 

= ^Yj «"*" [al u (0 + ai,_ k (o] ■ (i68b) 

k 

At the lowest-order approximation, by using the Bogoliubov transformations 

ao,k = «o,kVk + vo,kBJ_ k , (169a) 
«i,k = wi,kKk + v -i,k^li ,-k' (169b) 
«-i,k = «-i,k&-i,k + VijkSj _ k , (169c) 

the effective Hamiltonian for the noncondensate part can be written as 

*# = 2 ha>o^b ] 0k b . k + 2^ (fuoi,kb' l k t>i,k + fi^-i.kb[ lk b-i,k) , (170) 

k*0 k 

where £>o;±i,k are the annihilation operators of quasiparticles. As seen in Eq. ( I64ab . there is a two- 
fold degeneracy in energy hcji k = 7*&>-i,k due to the symmetry between the j — +1 sublevels. 
The ground state is defined as the vacuum of annihilation operators »o;±i,k' 

So,klg> = 0, (171) 
6±i,klg> = 0. (172) 
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From Eqs. (1168b and ( 1169b . the transverse magnetization density operators F+(r, t) and Fl(r, f) 
can be written in terms of quasiparticle operators £>o;±i,k as 

F+(r, t) = V2n^ ^ e~' k ' r [("* )k &I jk + v -i,k^-i -k) + ("-1 -k&~i ,-k + v i.-k&i jk )] 
k 

= V^JV*' [(*** + Vl *H,k + ( V -Lk + u -i-k)b-i,- k ] , (173a) 

k 

P-(r, t) = V2«o ^V*" [( M -i.k^-i,k + v i,kk-k) + ("i ,-k^i ,-k + v_i k )] 

k 

= A ^Je" a " r [(a* j k + v_i jU ) Slj k + (v* 1M + u lM ) hi,. k ] . (173b) 

k 

Here, we used M ± i jk = u ± \- k , v±i,k = v±i _ k for coefficients of the Bogoliubov transformations. 

If a particle with momentum fiko and spin j — 1 is created above the ground state, the system 
is in an excited state given by 

l e > =4kjg> = + "Ii^-L-H,] lg> = «I*$jJg>- (174) 

Using Eqs. ( 1173b and ( 1174b . we can calculate the equal-time spatial correlation of transverse 
magnetization with respect to this excited state as follows: 

(e\F + (r,t)F-(r',t)\e) 
=2« J] e" , ' (k ' r+k '' r °<e| [(< k + v lM ) b\ k (t) + (v*_ l k + M _ 1)k ) &_j _ k (f)] 

k,k' 

x [( M -i,k< + v-i,k')b f _ lk Xt) + (v* 1Js; + «i,k')^i,-k'(f)] |e) 
=2« J] e -' (kr+k ' r '>[ («^ + v ljk ) (v^ + « lik .) «^^^|«i*,| 2 W-**> 

k.k' L 



+ (v! 1(k + K_ ljk )( M ! 1(k , + v_i, k -)e- ! ' (<J - I - | '- w - 1 *' )f <5_ k j £ - 
=2« K, ko + vi, ko | 2 |Mi, ko | 2 ^ , ' ko ' (r_r ' ) + 2n J] l v -i,k + "-Lk| 2 e_ ' k 



-/k-(r-r') 



(175) 



Here, the last term in Eq. (1175b also exists in the spatial correlation of transverse magnetization 
for the ground state. Similarly, we have 



<eF_(r, t)F + (r', t)|e> =2« |< ko + v lM f \u lMo \ 2 e ik ^- r,) 
+ 2no2h'*i.k + "-i,k| 2 e ,1 " (r - r ' ) . 



(176) 



The equal-time spatial correlation of transverse magnetization with respect to the plane-wave 
excited state |e) then takes the following form: 

<e|F ± (r, t) ■ F ± (r', f)|e> - <g|F ± (r, t) • F x (r', f)lg> 

1 <eF + (r, f)P.(r', t) + fi_(r, t)F + (r', r)|e> - <g|F + (r, t)F_(r', t) + P.(r, t)F + (r', t)\g) 



=2n \u* l ko + vi iko | 2 |Mi, k „| 2 cos [k ■ (r - r')] . 



(177) 
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From Eq. ( 1177b , it is clear that, similar to the elementary excitation given by Eq. ( 1153b for the fer- 
romagnetic phase, the elementary excitation given by Eq. (1174b for the polar phase also shows a 
periodic spatial modulation in the system's transverse magnetization, and thus, can be interpreted 
as a transverse spin wave. We can, therefore, produce a localized wave packet of transverse mag- 
netization by locally exciting atoms initially in the j — to the j = 1 state (or equivalently, to the 
j = — 1 state). Due to the nonlinear dispersion relation [see Eqs. ( 1 1 39ab and ( 1140b l. the prepared 
spinor wave packet expands in space during its propagation. By measuring the group velocity or 
the rate of expansion of the wave packet, we can obtain information about the effective mass of 
the corresponding magnons. In the case of polar phase, the low-momentum region for which the 
dispersion relation of magnons is in a quadratic form is given by ej,' <K |ci|«o [see, for example, 
Eq. d 1 39ab l . The width of the initially prepared spinor wave packet, therefore, should be of the 
order of 10 //m, which can be produced by using a pair of laser beams whose beam waist is larger 
than that used in the ferromagnetic phase. Furthermore, in contrast to the ferromagnetic phase, 
the effective mass of magnons for the polar phase, and in turn, the time-dependent width of the 
spinor wave packet depends on the external magnetic field via the quadratic Zeeman effect as 
given by Eq. ( 1143b . A small variation of qs near 2|ci|«, which corresponds to a magnetic field 
of the order of hundreds milliGauss, can make a big difference in the dynamics of a spin wave, 
and thus, the time evolution of the spinor wave packet can also be exploited to perform a precise 
measurement on a magnetic field. 

6. Conclusion 

In this paper, we have studied the effect of quantum depletion at absolute zero on the energy 
spectra of elementary excitations for a Bose-Einstein condensate (BEC) of 87 Rb atoms in the 
F — 1 hyperfine spin manifold. We have generalized the Beliaev theory, which is a diagram- 
matic Green's function approach, to describe a system with internal degrees of freedom. The 
investigation was done on an atomic system whose ground state is in one of the two charac- 
teristic quantum phases of an F — 1 spinor BEC: the fully polarized ferromagnetic phase and 
the unmagnetized polar phase. In contrast to a spinless BEC, there are spin-wave elementary 
excitations in a spinor BEC in addition to the conventional density-wave excitation. We showed 
that the corresponding magnons in a spinor BEC have quadratic dispersion relations as opposed 
to the linear dispersion relation of phonons. We also found that in both cases of ferromagnetic 
and polar phases, the quantum depletion leads to an increase in the effective mass of magnons, 
while it does not alter the energy gap at the leading order. Under an external magnetic field, the 
effective mass of magnons for the polar phase depends on the strength of the quadratic Zeeman 
energy relative to the spin-exchange interatomic interaction, as opposed to the ferromagnetic 
phase. This demonstrates the difference in the coupling of the motion of magnons to the external 
field for the ferromagnetic and polar phases. Nevertheless, we found that the enhancement factor 
of the effective mass of magnons due to the quantum depletion turns out to be the same for both 
phases, and also independent of the external parameters of the system. This implies a universal 
mechanism whereby the quantum depletion affects the motion of magnons in spinor Bose gases: 
the motion of magnons is hindered by the interaction with the quantum depleted atoms. Further- 
more, in the case of 87 Rb atoms where the spin-conserving interaction is much larger than the 
spin-exchange one, the lifetime of magnons becomes much larger than that of phonons. This 
agrees with the mechanism of Beliaev damping as due to collisions between quasiparticles and 
the condensate, and can be understood by considering the momentum and energy conservations 
in the collisions. 
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For a system of ultracold atoms with a particle density n ~ 10 15 cm 3 , the effective mass 
of magnons increases by a factor of 1/[1 - 49 V«a 3 /(45 yjn)] ~ 1.01. The increase is about 
1%, which is expected to be measurable with high-resolution experiments. Moreover, by using 
a technique to effectively increase the scattering length a of the atoms, the quantum effect can 
become much larger, and easily measurable. To measure the effective mass of magnons in a 
dilute ultracold spinor Bose gas, we have proposed an experimental scheme which exploits the 
effect of a nonlinear dispersion relation on the spatial expansion of a spinor wave packet during 
its time evolution. By measuring either the group velocity or the rate of expansion of the wave 
packet of transverse magnetization, the information about the magnons' effective mass can be 
obtained, from which the quantum depletion effect can be probed. We also evaluated the time 
needed for the spinor wave packet to expand to the entire atomic cloud, and it is well within 
the lifetime of BECs in experiments of ultracold atoms. Using the fact that the effective mass 
of magnons for the polar phase is a function of the magnitude of external magnetic field, this 
kind of measurement can be used for numerous practical applications as, for example, to identify 
spinor quantum phases, or to be used for precision magnetometry. 
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Appendix A. Relation between T-matrix and vacuum scattering amplitude 

The T-matrix T^{p\,p2\pj L j>\) in the spin channel T, which is given by Eq. (l42l . satisfies 



the Bethe-Salpeter equation 13811 : 

hj (27!-) 



/d 
— |j V r (q)G ( Pl - q)G°(p 2 + q) 

xT r (pi-q,p 2 + q;p3,p4)- (A.l) 



This equation is illustrated in Fig. IA. 131 Because the form of the Bethe-Salpeter equation is the 
same for two spin channels T — and 2, the subscript T will be omitted below. 

Let us introduce the four-vector total momentum HP = hp\ + hp2 = hpj, + hp\, where the 
second equality indicates the conservation of total momentum and energy, and the four-vector 
relative momentum hp = (l/2)(ftpi - hp2), hp' = {l/2)(hp^ - Tip<\) for a pair of scattering 
particles. Equation ( 1A.U can then be rewritten as 

T{p, p', P) =V(p - V >) + if^f^SL V(q)G°(P/2 + p-q) 

x G°(P/2 -p + q)T(p - q, p', P), (A.2) 
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Figure A. 13: Bethe-Salpeter equation for the T-matrices Tj^(pi,p2; P3, p\) in spin channels T = and 2 [see Eq. jA.lH . 
The squares represent the T-matrices, while the free propagators, which describe spinless non-interacting Green's func- 
tions G°(p), are represented by solid lines with arrows. The wavy lines show the interatomic interactions Vy(p) in spin 
channels T = and 2. 



or in a form of an infinite series as 
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V(p, p', P) =V(p - p') + l - j £± V(q)V(p - q - p ) 

x j + p - q , P/2 + p - q)G°(P /2 - p + qo, P/2 - p + q) 

+ ■ ■ ■ . (A.3) 



Via a transformation of variables: go = qo + Po, the integral inside the square brackets in Eq. (1A.31 
can be rewritten as 

j |^GVo/2 - qo, P/2 + p - q)GVo/2 + qo, P/2 - p + q), (A.4) 

which is independent of po. In a similar manner, the higher-order terms represented by the dots 
in Eq. (lA.3t are shown to be independent of po and p' by iteration. Therefore, the T-matrices are 
independent of po and p' Q , and can be written as T(p, p', P). 

Next, we introduce a quantity ^(p, p', P) as an integration kernel of T(p, p', P) lfl H3~5ll : 

r(p,p',P)= f i!^V(q)^(p-q,p',P). (A.5) 

Note that Eq. ( IA.5t is similar in form to the equation relating the vacuum scattering amplitude 
-Mf(k, k')/(4-nh 2 ) to the scattering wavefunction iAk(p) in momentum space: 

— ^F(q)^ k (k'-q). (A.6) 
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From the Bethe-Salpeter equation [Eq. (IA.2H for F(p, p', P), we obtain the equation for;^(p, p', P) 

as 



Xip, p', P) =(2;r) 3 (5(p - p') + l - J ^G°(P/2 + p)G°(P/2 - p) 

3 V(q) X (p-q,p',P). (A.7) 



d 3 q 



(2tt) 

Indeed, by substituting Eq. ( IA.7t into Eq. ( IA.5I ). it can be seen that Eq. (IA.2b is satisfied. Cal 

p 



culating the integral with respect to po in Eq. ( IA.7I ) by using G a (p) = l/(po - + p + irj), we 



obtain 
X(P 



,p',P) = (2n) 3 6(p-p') + — — f A±V(q) X (p-q,p',P). (A.8) 



Note that Eq. (IA.8b for^(p, p', P) is similar in form to the Schrodinger equation for the scattering 
wave function </<k(p) in momentum space: 

^(p) = (2^(p-k)- 1 f ||v(q)ifk(p-i|). (A.9) 

Then, by using Eqs. (IA.6I ), dA.8b and dA.9l ), ;^(p, p', P) can be expressed in terms of i^k(p) and 
/(k', k) as (see, for example, 13511 ) 

xM , P) . Mp) + J ii fr(p)( ____L____ 

1 



tfq 2 _ ££i _ 

M M "I 



)/(p',q)*. (A. 10) 



Substituting Eq. ( 1A.10I ) into Eq. dA.5b . we obtain the expression of the T- matrix F(p, p', P) writ- 
ten in terms of /(k, k') as follows: 



r d 3 

r(p )P ', J p) = /(p,p') + J ^/(p.q)( 



i 



HPo-^+lp-Z^ + ir, 
1 



fr 2 q 2 fr 2 p' 2 

~M M ir l 



)/(p',q)*. (A.ll) 



From Eq. ( IA. lib , we can see that the T-matrix r^(pi, p2\ pz,p<C) = T^(p, p', P) can be fully ex- 
pressed in terms of the vacuum scattering amplitude -Mf^{p, p')/(47rh 2 ) in spin channel T ■ This 
scattering amplitude is a well-defined physical quantity even for a singular interaction potential. 

In the discussion of the T-matrix r j m j' m >(pu pi; Pz,Pa) in Sec. 12.31 we have neglected the 
dependence on the spin of intermediate states via the quadratic Zeeman energy qsij" + m") in 
the denominator of Eq. (l39l l and we are now in position to justify the validity of that approxima- 
tion. From Eq. d39l . the difference of Tj m ,fm'{p\,P2\ P3,Pa) between the cases in which the term 
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qsij" + m ") i s an d i s not neglected has the following order of magnitude: 

d 3 q qB 

(2^) 3 (2 M - 2e° q + ir,) fa - 2e\ - q B + irj) 

~ \c\\n V«a 3 

-CcomVho 3 (A. 12) 

Here, we consider only atoms in the low-momentum region er , , ^ , e° 4 <K co«, subject to a 
small external magnetic field q B ~ \c\\n «; cqti as discussed in Sec. 12.31 We also used /i ~ cqh, 
VV(p = 0) ~ co. From Eq. (1A.12I I. it can be seen that up to the order of magnitude of cqh V«a 3 , 
the approximation used to evaluate the T-matrix Yj m ,j'm'{p\,P2\ Pi, Pa) in Sec. I2.3l is justified. 



V r (v = 0) 



J 



Appendix B. Derivation of Eq. (l69l 

The second-order contribution to the T-matrix F^(p, p', P) given by Eq. ( 1431 is calculated to 



be 



r^(p, p', P) = Im{/ r (p, p')} + fr ( ^ 



h 2 n 2 h 2 p' 2 

~M M lT l 



AM ' M 

(B.l) 



where we neglected the momentum dependence of the generalized vacuum scattering ampli- 
tude /jr(p, p') in the q-integral, and replaced them with their zero-momentum limit fy-. These 
replacements are justified by the fact that the q-integral in the T-matrix contains f^, and is a 
second-order correction. 

From Eqs. (l53l l and (f4TT >. it can be seen that the contributions to the self-energies and chem- 
ical potential from the first-order diagrams in Fig. |5]involve the T-matrices F^(p/2, ±p/2,p), 
r^r(p, 0, 0) = Fyr(0, p, 0), and r^(0, 0, 0), whose second-order contributions are given by using 
Eq. jB.lt as 



f (p/2, ± p/2,p)=Im { /r(P/2. ± p/2)} + ^ f 0( — — 1— ^ 



h 2 t{ 2 h 2 p 2 

~~M AM lT l 



(B.2a) 



AM I* M 

1 



F^(p, 0, 0) = r£>(0, p, 0) = f r f |SL ( J q2 ■ + ^1 (B.2b) 

Here, we used the fact that the imaginary parts of the on-shell scattering amplitudes /^(p, p') with 
|p| = |p'| give second-order corrections [see Eq. (1711 1. while the off-shell scattering amplitudes 
fr(p, 0) and fr(0, p) are real numbers yj. 
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The q-integral in Eq. (lB.2ab can be rewritten in a form that is useful for the calculations in 
Sec. |4]by making a transformation of variables q = q' + p/2, with which we have 



k= q-p =q' -p/2, (B.3a) 

e o + e o = ^V + ^l, (B . 3 b) 
q M AM 

fi 2 n 2 fiV 2 

e o _ e o _ e o = _ (B 3c) 

p q AM M y ' 

jd 3 q = j d 3 q'. (B.3d) 



The q-integral in Eq. (IB.2al > then gives 
d 3 q 



J (2nf\ 



?- h 2 p 2 , ~ h 2 q 2 , . h 2 a 2 h 2 p 2 

foPO 4M M M - 1M - "1 



ap + 2yK - - + it] e$ - e% - + z?7 



(B.4) 



By substituting the lowest-order chemical potential in Eq. J53ct into Eqs. (1B.21 and (B.4\ . and 
using Eqs. d53l and flTt . we obtain Eq. ( 1691 . Here, for the second-order correction under con- 
sideration, the spin-exchange interaction ci(p, p') = [fz(p, p') _ /o(p, p')]/3 is neglected because 
of its small contribution compared with the spin-conserving one co(p, p')- 



Appendix C. Contributions to the self-energies and chemical potential from the second- 
order diagrams 

Appendix C.l. Ferrromagnetic phase 

In the second-order contributions to the self-energies and chemical potential, the spin-exchange 
interaction is neglected since it is smaller than the spin-conserving one by a factor of two hun- 
dreds. Therefore, the T-matrices in the second-order diagrams in Figs. [6]|9] are reduced to 

^jmj'm' — jj'&mm' i (C.l) 

where cq is given by Eq. (l49l i. On the other hand, the propagators in Figs. |6]|9] which are used 
to evaluate the second-order self-energies and chemical potential, are given by the first-order 
Green's functions in Eq. ( T571 i. Then, the contributions to the self-energy j}}.,{p) from the second- 
order diagrams in Fig. [6] are given as follows: 

( al ) = ^2"° C Yj ) W^A G >nm(q)G mm (q - p)6j t lS ft \ 

^ 3 q / Aj.qB^t A lM B lq 



_ «oCq r d J ( 
_ «oCq C d 3 < 



) 3 \p - U)i A - Wi, k + IT) po + tOl, q + Wi jk - it] 
q / {^l,q,#l,k} {^l,k,#l,q 



) 3 \po ~ Wl,q - Wl,k + IT] PO + 0)l,q + Wl,k ~ »7 



Sj,i6 f ,u (C2) 
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where k = q-p and |A;, q , B j;k ] = A jA B jiL + Aj M Bj, q . Here, Go(q) G (q- p) and G-\(q)G-\(q- p) 
give zero contributions to the go-integral in the first line of Eq. (IC.2K and in deriving the last 
line we used the fact that the value of the integral in the second line does not change under the 
exchange of variables q and k. Next, 



i 2 r d 4 q 
( a2 ) = ^2«oc J -^-^Gn{q)Gn{q- p)Sj i \6 fi 



= (al), (C.3) 
(a3)=(a2)= (al), (C.4) 



i 2 f d\ 
(a4) = ^-noco J G n (q)G n (q - p)5j, 



h 2'^o J (2tt) 4 JJ M 
q £i,k 

) 3 po - W ,q - <L>l,k + IT] 

B\,k 

) 3 po - W-l,q - Wi,k + iT] 



:(al)+ — J (2^ — 



TT TT^T — ^Oy,-H>/,-i> (C5) 



(bl) = ^n cg J ^Glfa^te-p)^,, 

f d 3 q / 1 1 \ 
tz— -t Ci, q Ci,k ; - (C.6) 



(b2)=(bl), (C.7) 
(b3)=(bl), (C.8) 
(b4)=(bl), (C.9) 



(cl) = ^2 n ° c o J (2^)4 Gll ^) G n^ ~ &> 6 J.i 6 f,i 

n Q cl r d 3 q / (Ai, q ,Ci, k ) (jW^k) \ 

= 7TTT ~ + ~ \Oj,\6y,i, (CIO) 

2ft 2 J (2^) j \ po - OJ\ A - wi, k + in Po + wi,q + wi,k - "7/ 



(c2)=(cl), (C.ll) 
(c3)=(cl), (C.12) 
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i C* d^ 

(c4) = — n 4 J Gjj(q)G\](q - p)5j 



H , u "» c o f d 3 q (-d, t ) 

h z J (2nY p - w ,q - wi.k + Z77 



^ r d3 q 

2 J (2jt)3 



noc^ f d j q (-Ci, k ) 



S^iSf-u (C.13) 



/j4 



= (cl), (C.14) 



(d2)=(dl)= (cl), (C.15) 
(d3)=(dl)= (cl), (C.16) 



C* d 

(d4) = ^«„c 2 J _^G iV ( ? )G 2 ;^-/^ 7 

= (c4), (C.17) 



' /"* d^ 

(el) = ^n c 2 J — Jj [G^( 9 X?ii(p " ?) " G%q)G\{p - q)]s jf 



_ n c z r d J q 
"~J (2^) 



) 3 \S(po - «l,q - Wi, k ) + /?7 ^(PO + 0)\ A + 0»i, k ) 



"oCo f d J ( 



o _ o 

fc k 

q / A i,k 



%0 - fq -%. + 2(C + Ci)« + IT] 
,2 r a! 



) 3 I h (p - a>o, q - wi,k) + ii/ 
1 



"oCq f d J ( 



o ,o 
" e k 

q / A i,k 



.2 ^ ^3 



) 3 \ft(po - w_i, q - w 1;k ) + irj 

1 



<H i"V. i. (C.18) 



%o - e q e k + 2(c + Ci)« + z'?7 

Here, we should subtract a term containing non-interacting Green's functions given by Eq. 
from the contribution of diagram (el) to avoid double counting of the contribution that has al- 
ready been taken into account by the definition of the T-matrix and first-order diagrams in Fig. [5] 
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Similarly for the diagram (e2), we have 



/ C d 

(e2) = ^n<,cg J — L [G n (?)G n (p - q) - G?($)G?(p - q)]Sj,iS f , 
" 3 q / Ai, q Ai, k Bi q B hk 



noCQ r d J { 

~~~J 



1 



frpO - e q - e k + 2 ( C + Cl)«0 + »7 



• (C19) 



Next, 



en f d 3 q 

¥j (C20) 



where we have introduced the convergence factor e m ° with 77 — > +0, which results from the 
normal order of field operators in physical observables. Similarly, we have 

=I/H^^- (C21) 

By summing up Eq. d69ab and Eqs. dC2l- dC2Ti we obtain Eqs. d7?1>-f75T> for Ejl (2) (p). 

Next, the second-order contributions to the self-energy 2)j»(/>) from the second-order dia- 
grams in Fig.Qare given as follows: 



_noc 2 Q r d 3 ( 
"~J (2* 

_ «qCq r d 3 < 
~ 2ft2"J (2tt 



3 q / Ai >q Bi ik Ai >k Bi, q 



) 3 \p - <^l,q - Wi, k + J77 po + Wi, q + £wi, k - J>7 



3 q / {Ai, q ,Bi, k } {Al,k,Bl,qJ 



tf/ltf/r.1, (C.22) 



( a2 ) = ^2 n ° c o J (2^4 Gll( -^ Gl1 ^ ~P) S i' lS rt 

= (al), (C.23) 

(a3) =(a2) = (al), (C.24) 



(a4)=(a2)= (al), 
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(C.25) 



(bl) = ^n c 2 J G?M)G 2 Mq ~ p)Sj,iS f , 



3 Ci : qCi,k 



1 



Po - OJl.q - Wi, k + it] po + COl,q + wi, k - it] 



(b2)=(bl), 
(b3)=(bl), 
(b4)=(bl), 



(cl) = 



{^l,q,Cl,k} 



{#l,q> Ci,k} 



PO - «l,q - 0»i, k + *>7 PO + Wi, q + 0»l jk - 177 



ft 



(c2)=(cl), 



(c3) 



«OCq f d 3 q / (Bl > k,C'l,q} {Ai,k,Ci, q } 

2^J (2^> 



PO - Wi, q - Oli,k + it] PO + <Wi, q + 0<l jk - Z?7j 



(c4)=(c3), 
(c5)=(cl), 
(c6)=(cl), 
(c7)=(c3), 
(c8)=(c3), 



(dl) 



(2;r) 4 



p)G}?( 9 )^,i^, 1 



= (bl), 
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(d2)=(dl)= (bl), 



(C.39) 



i r d\ 
(e) ^ C °J (2n 
= £o f d 3 q 
h J (2;r) 3 



G*{ q y™- C -^G\{q)G\{- q ) 



-Cl.q + 



C0«0 



2e° - 2(c + ci)n - »7 



(C.40) 



Here, we should subtract a term containing non-interacting Green's functions given by Eq. 
from the contribution of diagram (e) to avoid double counting of the contribution that has already 
been taken into account by the definition of the T-matrix and first-order diagrams in Fig. [5] We 
also have introduced the convergence factor e ma with 77 — > +0, which results from the normal 
order of field operators in physical observables. By summing up Eq. d69b| i and Eqs. ( IC.221 >- 
( IC401 . we obtain Eq. CTJi for Ejf 2) (p). 

It can be shown by changing the direction of momentum from p to — p that the contributions 
to Y?\{p) from the second-order diagrams in Fig. |8]are equal to Eqs. (IC .22t -( IC.401 >. In fact, it can 



be shown that - EjiOO to all orders (see, for example, J35J). Finally, the contributions to 

the chemical potential p from the second-order diagrams in Fig. |9]are given as follows: 



(al) 



r d 3 ( 

J (2* 



G mm (q)e" 



£0 

n 



3 B Uq , 



(C41) 



(a2) 



= (al), 



J (2n) 



G n (q)e" 



(C.42) 



(b) j [G\l(q) e ™ - C J^G\{q)G»{-q) 

r d 3 ( 

J (2w: 



£0 



TV 



c «o 



2e" - 2(c + ci)n a - it] 



(C.43) 



By summing up Eq. d69ct and Eqs. ( 1C.41I >-( |C.43| |. we obtain Eq. d78l l for p 



a, 
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Appendix C.2. Polar phase 

In a manner similar to the case of ferromagnetic phase, the contributions to the self-energy 
E,l(p) from the second-order diagrams in Fig. |6]are given as follows: 



riQcl 



r d 3 ( 

J (2n 



Ai,qBl,k 



Al,k#l.q 



Ao.qfio.k 



^0,k#0,q 



PO - W , q - (JL> ,k + ^ PO + W ,q + Wo,k ~ V) 



_ «qCq r d 3 ( 
~~ J (2* 



(^l,q.« 



l.k 



PO - Wi <q - Wi, k + /77 Po + Wi, q + Wi ik - !77 



1 

+ - 



{Ao,q, -Bo.l 



{Ao,k, #0,q} 



2\p - W ,q - W(),k + '?7 PO + W(), q + O» ,k - "7 

i , r d 4 g 
(a2) = ^2 "oc J G 00 (q)Goo(q - p)S jfi 5 /j0 

_n cl r d 3 < 



Ao, q , Bo,k! 



(Ao,k, #0,q 



q_ 

) 3 \po - <W ,q - w 0,k + 2»7 PO + W ,q + W ,k ~ "7 



(a3)=(a2), 



(a4) 



A/,qB(),k 



«oCq f d J q 



A , k B 



Po - oj jA - cjqm + in Po + oj j, q + w ,k - in 



Sit 



(bl) = ^n cl J ^| i [Gg( g )GS(9-p) + Glf_ 1 (?X??, 1 _i(?-/») 
+ Gi 2 u (ty)G! 1 u (< ? -/ 5 )]^ / ,o 



+ 2Ci >q Ci : k 



f d 3 q 



Co,qC(),k 



1 



1 



1 



PO - W ,q - W(), k + IT] po + W0,q + W0,k ~ *?7 
1 

Oj,0Of,0, 



Po - w i>q - a>ix + irj po + wi, q + o»i, k - *'?7 



(b2) = ^2 "° C J" 4j* G 00(?) G 00(? - P) S jfl S j'$ 

_ "oCq r d 3 < 



3 Co.qCo.k 



1 



po - a>o,q - w , k + in Po + w , q + w ,k - in 
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(C.44) 



(C.45) 
(C.46) 



(C.47) 



(C.48) 



6j,o5 f ,o, (C.49) 



(b3) = (b2), 



(C.50) 



(M)=(b2), (C.51) 



i f 

(cl) = -^"oCq J ^-3 G Q0 (q)G l ^(q - p)8j t0 5j,fl 

3 q / {^0,q, Q),k} {-Bo.q, Co 



_"oCq r 



....... . . 

Ao<5/,o, (C.52) 

P() - 0> , q - WO,k + 111 P0+ WO.q + W ,k - IT] 



(c2)=(cl), (C.53) 
(c3)=(cl), (C.54) 



7T^4 G jJ^ G li(l ~ P^Jf 



i 2 r d 4 < 

= /?"° C °J (2* 
_ "oCg T d 3 < 
"~J (2* 



3 q / ^;,qCo,k B/qCo.k 



(C.55) 



) 3 \ - LOj.q - coo,k + it] po + to/, q + w ,k - it] ' 



/ 2 r 21 

(dl) = ^2"° c o J G oo(^)G 00 (^ - p)5 jfi 8 r ,Q 

= (cl), (C.56) 



(d2)=(dl)= (cl), (C.57) 
(d3) =(dl) = (cl), (C.58) 



(d4) = -^n c 2 J --^ Gjj(q)G 2 ^(q - p)6 jr 

= (c4), (C.59) 



( el) = ^2 n ° c o / [ G i>(?) G oo(p - q) - G°j(q)G° (p - q)]s jf 

_ »qCq f d 3 q / A ; - q Ao,k B /,q B o,k 

^ J (2;r) 3 



K(po-cli lq - w ,k) + in h [p a + w jA + w ,k) - »7 



ftP() - e q - e k + 2c 0"0 - ?Bi 2 + ^ 
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2 ) = n ° c o J 7^~vJ [Goo(q)Goo(p -q)- G°(q)G°(p - qj\d jfi 5 /j0 

_ »qCq f d 3 q / A , q A ,k fio.q^o.k 

h J (2tt) 



ft (pO - <^0,q - ^0,k) + ir} h (po + <W ,q + W 0>k ) - Z77 



ftpo - e q - e£ + 2c «o + i»7 



<5;,o<5/,o, 



(C.61) 



i tr^ C d q 



(C.62) 



(f2) 4c r — 

_c f d 3 q 

" h J (2^)3 ^ 



(C.63) 



By summing up Eq. d93ab and Eqs. dC44l i-( IC63l l, we obtain Eqs. d94j» and (|95]) for SjJ ( \p) and 
Zq ( ' (2) (/5), respectively. 

Next, the contributions to the self-energy l} 2 j,{p) from the second-order diagrams in Fig. |7] 
are given as follows: 



_ n c 2 r d 3 ( 
" ft 2 J (2^r 



Al.q^l.k 



Al.k^l.q 



PO ~ Ul,q ~ <^l,k + If] PO + OJl, q + COix + IT] 



0,q»0,k 



A(),kB(),q 



PO ~ W , q - U>0,k + "7 PO + W ,q + W ,k + in) 

2 - - 3 ~ r/ {A,, q ,Bi, k ) {Ai, k ,Bi.q} 



noc\ C d 3 ( 
~J (2* 



PO ~ W 1>q - + in PO + OJ hq + Wl,k + "7 

{Ao,q, So,k} {-^0,k> #0,q ; 



1 

2\p - « ,q - w 0,k + ^ PO + w 0,q + W ,k + ^ 



(C.64) 



i 2 r ^q 

(a2) = t^«oc — — - t Goo(tf)Goo(# - p)£/,o£/,o 
rr J (.ZTrr 

_ «oCq T d 3 q / f^o.q, fio.k 1 



q_ 

T) 3 \p - <^0,q - W ,k + /?7 



l(),k, #0,q 



PO + W ,q + W0,k + ^ 



6 j,0$f,0, 



(C.65) 



(a3) = (a2), 
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(C.66) 



(a4)=(a2), 



(C.67) 



C 

(bl) = ^nocg J ^[G^G^-^ + Gl^^^^-p) 
+ Gi 2 u (^)G 21 u ( 9 -/ 3 )]^ / , 

""^J (2a-; 



r) 3 



Co,qC(),k 



1 



+ 2Ci >q Cl i k 



PO - «0,q _ W 0,k + IT] PO+ W ,q + W ,k - IT]) 
1 1 



PO - OJl,q ~ U\,k + IT] PO + OJl A + OJix - IT] 



(C.68) 



(b2) = Jf n ° c o j G oo(?) G oo(? - P)8j,oSf,o 
_ n Q c 2 f d 3 q 



1 



PO - W ,q - U) ,k + irj Po + w , q + wo,k - "7 



^/,o, (C.69) 



(b3)=(b2), 



(C.70) 



/i4 
^ [ol?_ifo)G£fo - p) : + «>,. i5,.i) 

+ Gw(q)G 2 ^(q - p)6 jt oSffi] 



n c, 



J (2tt; 



Cl,qCo,k 



1 



1 



PO - <^l,q - W ,k + "7 PO + Wi, q + fc) ,k - ^ 
1 



PO ~ W ,q - O> ,k + IT] 



1 
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(C.71) 



(cl) = ^j"oc5 



2 f d 4 q 
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Goo(^)G 00 (^ - p)6jflSf,o 
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Bq a , Co,k\ 



) 3 \ po - W ,q - W ,k + "7 PO + W 0) q + W ,k - IT] 



StfSffl, (C.72) 
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i C d^ 

( c2 ) = ^2 n ° c o J [Goo(?)G! 2 _i(^ - p)(6 jt \S r -\ + 8j-\6j>,\ 

+ Goo(q)G l w(q - p)S jt0 6ffl] 
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(C.73) 
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+ Gooiq - P)G^(q)Sj t0 6ffl] 
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1 {fio,k,C ,q} A ,k,C ,q} 
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(C.75) 



(c5)=(cl), 



(c6)=(cl), 



(c7)=(c3), 



(c8)=(c3), 



(C.76) 
(C.77) 
(C.78) 
(C.79) 



f d^ 

(dl) = ^n cl J ^ [g£(* - p)G\ 2 _ l (q)(S M 6j^ l + Sj, „>>. : 



+ G oo(<? _ P)Gw(q)6jfldf t a 
= Cb4), 
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(C.80) 



(d2) = ^2 n ° c o J (2^)4 G oo^) G oo(? - P)$j,o 6 f,o 
= (b2), 



(C.81) 
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+ [G 5( ? ) e ' Wo - con G (q)G°(-q)]6j, S rfi } 
co f d 3 q [ 

/ c n 

+ " C ,q + 



2e° - 2c «o - it] 



(C.82) 



By summing up Eq. d93bb and Eqs. ( tC.641 h fC.821 l. we obtain Eqs. 
respectively. 



and (|57) for s! 2 ^ (p) and 



-oo 



As in the case of ferromagnetic phase, it can be shown that 1?^,{p) - ^)i(p) by changing 
the direction of the momentum from p to —p and using the spin symmetry for the polar phase. 
Finally, the contributions to the chemical potential fi from the second-order diagrams in Fig. 
are given as follows: 



n 



(C.83) 
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i f d\ 
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ft J (2tt 
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(C.84) 



C d^ 

(b)-^co J -^[G^(q)e^-c n G (q)G (-q)] 

= £ o r_d^ 

ft J (2tt 



^(-Co, q+ 20 C 2 °"° . )■ (C.85) 
By summing up Eq. d93cT > and Eqs. (03t - dC85l . we obtain Eq. (gSJ for fi (2) . 
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Appendix D. Imaginary parts of self-energies 



Appendix D.l. Ferrromagnetic phase: l} Q ^ 2 \p) 

By making a transformation of variables q = p/2 + q', we have 

k= q-p=q'-p/2, (D.l) 

4-4-4=2{4i2-4'), (d-2) 

J (2*)3 J (27T)3- ^ ^ 

The imaginary part of the last term in the second line of Eq. (l75l l can then be rewritten as 

( fo+Vj ) f d 3 q 1 
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^o 2 + 



2-| 



(D.4) 



This cancels with the first term in Eq. (175V Therefore, by limiting our consideration to a small 
external magnetic field qg ~ \c\\n <k cqh, and ignoring any difference of the order smaller than 

11(2), 



co« yna 3 , the imaginary part of E. |f) (p) is reduced to 



t v ii ( 2), , "oc / r d 3 q Ai, k + fli, 

ImX oo V) = ir Im 7Tl3 
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(D.5) 



We then have 
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d(cos 0) 



^i,k 



(5(w ,p - W(),p + k - to\,k)- ( D - 6 ) 



62 



Here, 6 is the angle between p and k. The argument of the Dirac delta function is 
wo, p - a>o,p+k - w^k 



e°-e° 

fc P fc p+k 



fr|p||k|cos<9 hk 1 



M 



2M 



J$<fl cos 6 - e£ - ^e°[e° + 2(c + Cl )n ] 



= -y]4 + 2(c + ci)« + ^4 + 2 cos 



//L 



(D.7) 



For the low-momentum region under consideration ej,' <K co«o, the expression inside the square 
brackets of the last line in Eq. (ID. 7b is always positive for any value of 9 e (0, 7r). Therefore, the 
argument of the Dirac delta function vanishes only at = 0, and the value of the integral in the 
last line of Eq. dD.6b is, to within a multiplying factor, given by 



lim 



6 k 



1 



hd>[ k ^(^O.P ~^0,p+k -<*) 1 .k ) 
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ci)n ] 



2£°+2(c +C])n 



- v 'e°[e°+2(ro+ci)no] 



+ 1 + 



yJf$co$6 



=0. 

This implies that 

Similarly, we have 



dltr^ 2 \p) 



= 0. 



(D.8) 



(D.9) 
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de k — — d(cos0) — 5'(w , P - w ,p+k - wi ik ), 

J 2 V2tt 2 J Swi.k 

-1 

(D.10) 
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where 5'{x) is the first derivative of the Dirac delta function. Using the identity 

w(*)]y 



S'[f(x)l 



fix) 
[6(x - x )/f'(x )]' 

fix) 
6'(x - x Q ) 



f'(xo)f'(xy 

where xq is the zero point of function f(x), we have 
SlmZ'fV) 
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= 0. 



(D.12) 



Here, the multiplication factor outside the integrals in Eq. ( ID. 121 ) corresponds to f'(xo) in Eq. (ID. 1 U . 
From Eqs. flD.9t and (1D.12I I. we have 



lm^ (2) (p) = + 0[(p -wo,p) 2 l 



(D.13) 



Appendix D.2. Polar phase: S[[' 2 '(p) 

The imaginary part of xj [ <2) (p) is given by 

n cl r d 3 q (-ne£) 



11(9 , «oc r d q v — /lc w r 
lmE n k \p) = J — ^ — — [Ai, q ^(p - wi >q - w ,k) + B hq 6(p + wi, q + w .k)J ■ 



(D.14) 
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Iml!! (2) (p) 
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(D.15) 



Here, is the angle between p and k, and in deriving the last line of Eq. (ID. 1 5I > we used the fact 
that the argument of the Dirac delta function <5(wi jP +wi jP+k +Wo,k) is always positive. We consider 
only the low-momentum region e p <zz \c\\n$ and the external parameter region q B +7.c\n§ ~ \c\\hq. 
Fork such that |p+k| > |p|, we have e° +k > e p , wi, p+k > wi iP , and, in turn, wi, p -wi iP+k -wo, k < 0. 
In contrast, for |p + k| < |p|, we have e ( ^ +k < e p <g \ci\no and |k| ~ |p|, e k ~ e p «: \c\\no «: cqhq. 
The argument of the delta function in Eq. (ID. 15b is then reduced to 



^l,p _ Wi, p+k - W(),k 

V?b(?b + 2ci«o) 



1 



y/qB(.qB + 2cin ) 



h 

q B + C\tiQ 



2\qB qB + lciiiQ 
1 



el + 



p 



2\q B q B + 2cm j P +k 



+ 



p 

\c\\n 



^ fc P fc p+k' 



V?b(9b + 2ci« ) ^ 



<^0.k 



= - a e 
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e° + 2c «o 



W0,k 



(D.16) 



Because <K co«o> the expression in the square bracket of the last line of Eq. (ID. 16b is always 
positive for any value of 6 e (0, n). Therefore, the integral in the last line of Eq. (ID. 61 is, to 
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within a multiplication factor, given by 

t^O V "kink 



lim Je. A\„ + \,— 
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This implies that 
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Using the identity ( 1D.1 U . we have 
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dpo 



oc lim 



1 



X 



1 



(^) 3/2 



7^ A[ p+ k 
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From Eqs. (ID. 18b and (ID. 20b . we obtain 



imz;;%) = o+tf[(p -wi.p) 2 ] 
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Appendix E. Real parts of self-energies 



Appendix E.l. Ferrromagnetic phase: l}^ 2 \p) 

By limiting our consideration to a small external magnetic field qs ~ \c\\n <k cqii, and 
ignoring any difference of the order smaller than \cq\h V«a 3 , which is justified at the second- 
order approximation, the real part of Y^^ip) given by Eq. ( l75l l is reduced to 
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where f denotes the principle value, and 5 is defined by Eq. (|5TT i. Here, we used Aik 
2Ci, k = e°/(fiwi, k ) and 
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f ^ B - 



(c + ci)n M 



2 V2ft 
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3tt 2 ft 
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3/2 V2 
3 



«()fl 3 , 



where we have ignored terms that contain the factor |ci|/co <z. 1. 

First, we calculate the third line of Eq. (IE. 11 . Putting q = p/2 + q', we have 



k = q-p 
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The third line of Eq (IE. lb then becomes 
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By taking a transformation of variables |q| — > e^, and using the indefinite integration 
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together with the definition of the principle value 
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For the remaining term in Eq. (IE.U , its value in the low-momentum region ej] <K co« is 
obtained analytically by making Taylor expansions around po = wo, P and p = as described in 
Eqs. ( l82l and ( l85l l. The expansion coefficients are then calculated as follows: 
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Appendix E.2. Polar phase: ZjJ (p) 

Neglecting terms of the order smaller than co«o y/noa 3 , which is justified at the second-order 
approximation, the real part of X[] <2) (/?) is then reduced to 
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Here, we used 
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First, we consider the first term in Eq. (IE. 16) : 
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(E.19) 



Here, as moving from the second line to the third line in Eq. ( IE. 19b we used a transformation of 
variables [see Eq. (1E.3H . Furthermore, in the third line the main contributions to the first and the 
second integrals arise from e q ' ~ \c\\n and e° ~ co«, respectively, which results in the fact that the 
first integral is smaller than the second one by a factor of the order of Vki|/co <*c 1, and thus, 
the second integral was neglected. The integral in the next to the last line was directly calculated 
by using 
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(E.20) 



where x x = lirn v _>oo, a ± = a + 6(6 — > +0). 

For the remaining term in Eq. (IE. 16) , we can obtain an analytic result for the low-momentum 
region e'' «C \c\\n <k cqii by making Taylor expansions around po = o)\ >p and p = as described 
in Eqs. (1102b and ( 1105b . The coefficients of the expansions are then calculated as follows: 
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ReZ[[ (2) (p = wi ;P ) 
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Here, we used the fact that the main contribution to the integral in the second line of Eq. ( 1E.22I > 
comes from e° ~ co«o » ci«o, and we can approximate ftwi q e^', ftwi p= o - 0, Ai q 1, Z?i q 
0. This is because the integral converges at both the upper limit » co«o, and the lower limit 



e q — > 0. Next, we have 



5ReE|J (2) (p = wi, p ) 



<3w . 



= 0, 



p=0 



(E.23) 



d(coo,p) 2 



p=0 



q B + ci« 



71 



V«o5 3 , 

«0C0 



I 3tt 2 jq B (q B + 2cin ) 3607r2 J 
where we used the result of the following integral: 

13jc 2 + 1 lx 3 - V* 3 (2 + x) + 29 V^ 5 (2 + x) + 8 V* 7 (2 + *) 71 



dx ■ 



3(2 + x) 5 ' 2 (x + \Jx(2 + x)) 
Therefore, the real part of the self-energy £[[ (2, (po = wi, P ) can be written as 



90 V2 



, 5 « ( )C / — rr , 

37T 2 ft 

x -y/iofl ■ 

"oco 



1 



9b + c l"0 



71 



6^ 2 V?b(?b + 2cm ) 7207r2 J 



(E.24) 



(E.25) 



(E.26) 



72 



In a similar manner, the real part of the self-energy 2j[ (2) (-/7)|p 0=a)1 can be calculated by replac- 
ing ksq-p with k = q + p, and we obtain 



Res;; (2) (- P )i P0= „, p 



5 HqCo I — TT 

; 3^— V " 0fl + 

x V«ofl ■ 

"oco 



q B + ci« 



71 



6^r 2 ^Jq B (q B + 2 Cl «o) 7 207r2 



(E.27) 



Next, we have 



5Rez;; (2) ( P ) 



Po=<<»i,p.P=0 



»qCq f d 3 q ^ 
H 2 J (2n) 3 ha>o, q 



(A>l,p=0 - Wi jq - Wo,q) 2 



B 



l.q 



(E.28) 

(fc>l,p=0 + Wi >q + W ,q) 

As above, the main contribution to the integral in Eq. (IE. 28) arises from e° ~ cqhq » c\Hq. 
We can then approximate tko\ A e?, fayi )P =o - 0, Ai, q 1, Z?i, q 0, and the integral can be 
evaluated straightforwardly as 



3ReS;; (2) (p) 



dpo 



Po=^i, P .P=0 



« c 2 M 3/2 
tt 2 V2ft 4 



e2de° 



,2 M 3/2 



Similarly, we have 



V"oa 3 - 

J7P- 



3Re2;; (2) (-p) 



■(«oco) 



q q w ,q(eq + ha> 0A ) 2 

oo 

1/2 j yfxdx 



y/x{x + 2)[x + y/x(x + 2)] 2 



<3/?c 



Po=^i, P ,P=0 



For the derivatives with respect to wo, P , we obtain 



dco , p \ dp 



= 0, 



(E.29) 



(E.30) 



(E.31) 



p=0 



9 


fdReL\f\-p) 




\ dpa 



P0=^l,p 



p=0 



d 2 / 3ReZ;; (2) (p) 
<5(wo, P ) 2 \ dpo 



/>0=«Jl.p 



1 gB + cjno 



p=0 



3^ 2 + 2cmo) 607f2 

ft 2 



(«oco) 2 ' 



73 



(E.32) 



(E.33) 



d 2 /dRe£ 



11(2)/ 



d(coo,p) 2 \ dp 



p) 






f 1 , H \ 

1 q B + ci« 7 




PO=^l,p^ 


p=0 


l 3?f2 V^(4B + 2ci« ) 60;r2 J 



/r 



(«oco) 2 



(E.34) 



From the above results, we obtain the real parts of the self-energies S[J (2) (+p) as given in 
Eqs. 033 and (fTTl . 
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